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Abstract 

We establish a new criterion for the dynamical stability of black holes in D > 4 
spacetime dimensions in general relativity with respect to axisymmetric perturba- 
tions: Dynamical stability is equivalent to the positivity of the canonical energy, 
£, on a subspace, T, of linearized solutions that have vanishing linearized ADM 
mass, momentum, and angular momentum at infinity and satisfy certain gauge 
conditions at the horizon. This is shown by proving that — apart from pure gauge 
perturbations and perturbations towards other stationary black holes — £ is non- 
degenerate on T and that, for axisymmetric perturbations, £ has positive flux 
properties at both infinity and the horizon. We further show that £ is related 
to the second order variations of mass, angular momentum, and horizon area by 
£ = 5 2 M — J2a ^a3 2 Ja — -^S 2 A, thereby establishing a close connection between 
dynamical stability and thermodynamic stability. Thermodynamic instability of 
a family of black holes need not imply dynamical instability because the pertur- 
bations towards other members of the family will not, in general, have vanishing 
linearized ADM mass and/or angular momentum. However, we prove that for any 
black brane corresponding to a thermodynamically unstable black hole, sufficiently 
long wavelength perturbations can be found with £ < and vanishing linearized 
ADM quantities. Thus, all black branes corresponding to thermodynmically un- 
stable black holes are dynamically unstable, as conjectured by Gubser and Mitra. 
We also prove that positivity of £ on T is equivalent to the satisfaction of a "local 
Penrose inequality," thus showing that satisfaction of this local Penrose inequality 
is necessary and sufficient for dynamical stability. Although we restrict our consid- 
erations in this paper to vacuum general relativity, most of the results of this paper 
are derived using general Lagrangian and Hamiltonian methods and therefore can 
be straightforwardly generalized to allow for the presence of matter fields and/or to 
the case of an arbitrary diffeomorphism covariant gravitational action. 
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1 Introduction 

It is of considerable interest to determine the linear stability of black holes in D-dimensional 
general relativity. It is also of interest to determine the linear stability of the correspond- 
ing black branes in (D + p)- dimensions, i.e., spacetimes ^ x T p with metric of the form 

d^ +P = d4 + I>, 2 . (1) 

i=i 

where ds 2 D is a black hole metric. One can analyze this issue by writing out the linearized 
Einstein equation off of the black hole or black brane background spacetime. One way 
of establishing linear stability is to find a suitable positive definite conserved norm for 
perturbations. Linear instability can be established by finding a solution for which some 
gauge independent quantity has unbounded growth in time. However, even in the very 
simplest cases — such as the Schwarzschild black hole (331 ED] and the Schwarzschild black 
string [in] — it is highly nontrivial to carry out the decoupling of equations and the fixing 
of gauge needed to determine stability or instability directly from the equations of motion. 
It is particularly difficult to analyze stability when the background is stationary but not 
static. In the case of a Kerr black hole in 4-dimensions, the Teukolsky formalism [38] 
reduces the coupled system of linear perturbation equations to a single equation for a 
gauge invariant, complex variable, but there is no known formalism of comparable power 
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to decouple the equations in other cases, including the higher dimensional Myers-Perry 
analogs of Kerr; see [H El [7] and references therein. Thus, it would be extremely useful 
to have a stability criterion for black holes and black branes that does not require one to 
perform a complete analysis of the linearized perturbation equations. 

For the case of black branes, a simple criterion for stability was proposed in [TBI HH], 
based on the analogy between laws of black hole mechanics and the laws of thermody- 
namics. We now describe their criterion. Let us assume for simplicity that we are given 
a stationary, asymptotically flat black hole solution with compact event horizon cross 
section, in vacuum general relativity (no matter fields) in D spacetime dimensions. By 
the rigidity theorem [27J , we may assume that our solution is not only stationary and 
asymptotically flat, but has additional commuting rotational Killing fields 
with associated angular momenta Ji, . . . , Jn, where N < |_^y^J • Suppose now that we do 
not just have one such solution, but — as usually happens — a family of solutions parame- 
terized by the mass M, and the angular momentslj Ji, . . . , Jn, such as the Myers- Perry 
family [21]. As usual, we identify the entropy as S = \A, where A = A(M, Ji, . . . , JjA is 
the area of the horizon. Consider the Hessian matrix of S, 



Then by analogy with the thermodynamical properties of ordinary laboratory type sys- 
tems, the signs of the eigenvalues of Hesss might be expected to be related to the stability 
of the system: For normal laboratory-type systems we expect the system to be unsta- 
ble if the Hessian has a positive eigenvalue, which would correspond to a negative heat 
capacity. Therefore, one might expect that negativity of the Hessian, eq. (EJ, would 
be a necessary condition for black hole stability. However, this is clearly false: For ex- 
ample, the 4-dimensional Schwarzschild black hole has negative heat capacity, since its 
area is A = 16irM 2 , so d 2 S/dM 2 > 0. However, despite its negative heat capacity, the 
Schwarzschild black hole is well known to be linearly stable. Thus, for black holes, ther- 
modynamic stability with respect to perturbations to other stationary black holes cannot 
be necessary for dynamical stability. Nevertheless, there is no obvious counterexample 
to the possibility that thermodynamic stability could be necessary for black brane sta- 
bility, since, in particular, the Schwarzschild black string is known to be unstable [T6] . 
Therefore, suppose we exclude consideration of black holes and restrict attention to black 
brane spacetimes of the form ([TJ, where the Zi coordinates are periodically identified with 
length I, so that horizon cross-sections are compact and have a well defined area. The 
Gubser-Mitra conjecture states that any such black brane spacetime is unstable for suf- 
ficiently large / if the Hessian matrix (of the original D-dimensional black hole family) 
has a positive eigenvalue. In this paper, we will show that this conjecture follows as a 
consequence of a more fundamental dynamical stability criterion that we shall establish. 

Another simple possible stability criterion that is applicable to black holes is the "local 
Penrose inequality," discussed recently in [10J. We reformulate this criterion as follows: 
Suppose, as above, we have a family of stationary, axisymmetric black hole spacetimes 
with mass M and angular momenta Jj, . . . , J^. Since surface gravity cannot be nega- 
tive, it follows immediately from the first law of black hole mechanics that the horizon 

In the presence of matter fields, there may be additional conserved charges. 





3 



area, A, must be a non-decreasing function of M at fixed J\, . . . , Jn- Let 5W>(A) be an 
arbitrary smooth one-parameter family of axisymmetric metrics with g a &(0) being equal 
to a member of this family with non-zero surface gravity. Let B denote the bifurcation 
surface of the Killing horizon of g a b(0)- Let .A(A) denote the area of the apparent horizon 
of S'aft(A) on an initial slice passing through B. As will be shown in subsection 2.1 be- 
low, to first order in A the event horizon coincides with the apparent horizon. To second 
order in A, the event horizon need not coincide with the apparent horizon — assuming 
cosmic censorship, it cannot lie inside the apparent horizon but it may lie outside of it. 
Nevertheless, to second order, the only contribution to the difference between the area, 
A, of the event horizon and A arises from the (second order) displacement of the event 
horizon relative to the apparent horizon in the background metric g a b(0)- However, in the 
background metric, B is an extremal surface^). It follows that to second order, we have 
A(X) = A(X). Now, let g a &(A) be the one parameter family of metrics in the black hole 
family such that M(A), Ji(A), . . . , Jjv(A) agrees with the corresponding ADM mass and 
angular momenta M(A), Ji(A), . . . , Jat(A) of g a bW- Then, by the first law of black hole 
mechanics together with the fact that the bifurcation surface of g a b{fy has vanishing ex- 
trinsic curvature, we automatically have (dA/d\)\\ = o = (dA/d\)\\=o- However, suppose 
that S 2 A = (d 2 A/d\ 2 )\x=o > 5 2 A = (d 2 A/d\ 2 )\\=o- Then, since horizon area is nonde- 
creasing with time (assuming the validity of cosmic censorship), M is non-increasing with 
time (by positivity of Bondi flux), and J a is conserved (by axisymmetry) , we obtain a 
contradiction with the possibility that g a b(A) could "settle down" at late times to a black 
hole solution in the original family. This strongly suggests that satisfaction of the "local 
Penrose inequality" 

5 2 A<5 2 A (3) 

should be a necessary condition for stability of any black hole member of a family. This 
argument does not suggest that satisfaction of the local Penrose inequality should be 
sufficient for stability with respect to axisymmetric perturbations, but evidence in favor 
of this conjecture has been presented in [10]. In this paper, we will prove that our more 
fundamental stability criterion implies that satisfaction of the local Penrose inequality is 
necessary and sufficient for black hole stability with respect to axisymmetric perturba- 
tions. 

We now describe our main results. We consider stationary, asymptotically flat black 
holes in D > 4 spacetime dimensions in vacuum general relativity. In section 5, we will 
also consider the corresponding black branes in D + p dimensions. By the rigidity theo- 
rem [2~Tll2~T] . the event horizon of the black hole must be a Killing horizon, and, if rotating, 
the spacetime must be axisymmetric (possibly with respect to multiple rotational planes 
if D > 5). We restrict consideration to the case where the surface gravity, k, of the event 
horizon is non-vanishing. In the case of a nonrotating black hole, our results apply to 
arbitrary, smooth, asymptotically flat solutions of the linearized field equations, but for a 
rotating stationary, axisymmetric black hole, some of our results apply only to perturba- 
tions that preserve axisymmetry in the rotational planes that have nonvanishing angular 
velocity, Qa, of the horizon. For perturbations of the black hole, we define a quantity 
£ — called the canonical energy — which is quadratic in the perturbation and given by an 
integral over a Cauchy surface, E, for the exterior region. (Thus, £ extends from spatial 

2 This follows from the fact that both its outgoing and ingoing expansion vanish. 
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infinity to the bifurcation surface, B, of the black hole.) We will establish the following 
properties of £: (1) Without loss of generality, certain gauge conditions can be imposed 
near the horizon of the black hole that, in particular, ensure that the location of the true 
event horizon does not change to first order. When restricted to axisymmetric perturba- 
tions with vanishing linearized linear momenta, 5Pi = 0, and vanishing linearized change 
of horizon area, 5 A = 0, the quantity £ is invariant under axisymmetric gauge transfor- 
mations that satisfy these conditions at the horizon and approach arbitrary asymptotic 
symmetries at infinity. (2) £ is conserved in the sense that it takes the same value if 
evaluated on another Cauchy surface E' extending from spatial infinity to B. (3) For ax- 
isymmetric perturbations £ is related to the second order variations of mass, M, horizon 
area, A, and angular momentum, J A byEI 



(4) A Hilbert space V can be defined whose elements are square integrable initial data 
that (weakly) satisfy the linearized constraints, satisfy our horizon gauge conditions and, 
in addition, at spatial infinity have vanishing linearized ADM mass, linear momentum, 
and angular momentum with respect to the rotational Killing fields of the background 
spacetime. There is a dense domain, T, of smooth elements of V on which £ can be 
defined as a quadratic form. £ is symmetric onTxT and is degenerate precisely on the 
perturbations that are "perturbations toward stationary black holes." Thus, if we define 
a new space T' by modding out by perturbations toward stationary black holes (which, 
of course, includes all pure gauge perturbations), £ is non-degenerate. Consequently, on 
this space, either (a) £ is positive definite or (b) there is a \& e T 7 such that £('$, < 0. 

(5) Consider a foliation by Cauchy surfaces, E(t) for the exterior region such that each 
E(t) is composed of a hypersurface, y(t), that extends from a cross-section C{t) of future 
null infinity to a cross-section B(t) of the black hole horizon together with the 
portions of future null infinity and the future horizon that lie to the past of C(t) and 
B(t) respectively. For perturbations that preserve axisymmetry in the rotational planes 
that have nonvanishing angular velocity, Qa, of the horizon, the canonical energy on E(£) 
(drawn as the broken red line in the figure) can be expressed as the sum of a manifestly 
non- negative contribution from (equal to the Bondi energy flux), a manifestly non- 
negative contribution from (given by an integral involving the square of the linearized 
shear), and a term £{t) which differs from the canonical energy evaluated on y(t) only 
by boundary contributions from C(t) and B(t). Furthermore, the boundary contribution 
from C(t) vanishes if the Bondi news vanishes on C(t) and the boundary contribution 
from B(t) vanishes if the perturbed shear of the horizon vanishes on B(t). 

3 Although S 2 M, S 2 A, and 5 2 J a individually depend upon the second order metric perturbation, it 
follows from the first law of black hole mechanics that the combination appearing on the right side of (@| 
is independent of the choice of second order metric perturbation. 

4 In order to make use of the machinery of null infinity, we must restrict consideration to even di- 
mensional spacetimes, since J? does not exist in odd dimensions [23] ■ We do not believe that this is an 
essential restriction, i.e., we believe that our results on the positivity of flux at infinity will also hold in 
odd dimensions. 




(4) 
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A stability criterion for black holes and black branes with respect to axisymmetric 
perturbations can now be formulated as follows: Suppose that £ is positive definite on T' 
(case (a) of point (4) above). Since £ is conserved, we obtain a positive definite conserved 
norm on T' . Since T' excludes from T only those perturbations that are towards station- 
ary black holes — and stationary perturbations are automatically stable — this establishes 
stability of perturbations in T, at least in the sense of establishing the nonexistence of 
"growing modes". Furthermore, if the black hole or black brane whose stability is being 
investigated is a member of a familjQ of stationary black holes/branes containing arbi- 
trary variations of M, J\, . . . , J/y, then a general perturbation can be written as the sum 
of a perturbation in T and a perturbation to another member of the family, thereby 
establishing stability for general perturbations. 

On the other hand, suppose £ fails to be positive definite on T 7 and hence, according 
to point (4) above, there exists a perturbation in T' for which £ < (case (b)). We obtain 
a contradiction with the stability of the black hole or black brane by the following type 
of argument used previously in [HJ [121 H31 EH] : By property (5), the quantity £{t) is non- 
increasing with time, and, since £ = £ on S, for this perturbation we have £{t) < £ < 0. 
If the black hole/brane were stable, then the Bondi flux at J? + and the perturbed shear at 
Jf? + should go to zero at late times for this perturbation. Hence, £{t) should approach the 
canonical energy evaluated on S^{t). However, we have just seen that the canonical energy 
is bounded away from zero on y(t) at late times. Consequently, the perturbation cannot 
be approaching a solution that is a perturbation towards a stationary black hole since 
the linearized mass, linear momentum, and angular momentum must remain zercH and £ 
vanishes for perturbations towards a stationary black hole with vanishing linearized mass, 
linear momentum, and angular momentum. On the other hand, as already stated, any 
stable perturbation should become "non-radiating" at late times. However, this should be 
impossible, since non-stationary, non-radiating linearized perturbations should not exist. 

Our argument for stability in case (a) is a genuine proof that perturbations remain 
bounded in a suitable norm, although it does not establish pointwise boundedness and 

If the black hole or black brane in question is not a member of a general family, one would have to 
separately establish the existence of a stable perturbation for given variations of M, J\ , . . . , Jjv ■ 

6 Mass and momentum can be radiated at quadratic order; angular momentum cannot be radiated at 
all for axisymmetric spacetimes. 
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does not establish decajo On the other hand, our argument for instability in case (b) has 
a status closer to that of a plausibility argument than a complete and rigorous proof. At 
the end of subsection 4.2 below, we shall indicate some possible strategies for improving 
upon these results. 

At a conceptual level, our above criterion for black hole and black brane stability is 
remarkably simple: Evaluate the canonical energy £ of initial data for perturbations that 
(i) satisfy the constraints, (ii) have vanishing linearized ADM mass, momentum, and an- 
gular momentum at spatial infinity, and (hi) satisfy our gauge conditions at the horizon. 
If £ > for all such perturbations, then the black hole/brane is stable; if not, then it is 
unstable. Thus, £ provides a "variational principle" for black hole and black brane stabil- 
ity, and one can test for stability by plugging in "trial initial data" into the formula for £. 
Of course, the trial initial data must be chosen so as to satisfy the linearized constraints, 
so this variational principle cannot be used as readily as unconstrained variational prin- 
ciples^). Nevertheless, the use of our variational principle — even if it requires solving the 
constraints — represents an enormous simplification over performing a complete analysis 
of the dynamical behavior of solutions to the full set of perturbation equations. 

If we have a family of black holes or black branes parameterized by M, J±, . . . , Jjv, then 
we shall show at the beginning of section 5 that positivity of the right side of (jlj) for all 
second order variations within this family is equivalent to the negative definiteness of the 
Hessian (12"1). Thus, thermodynamic stability of the family as defined above by the Hessian 
criterion is equivalent to positivity of £ for perturbations within the family. However, 
our criterion for dynamical stability is positivity of £ for the case of perturbations with 
vanishing linearized ADM mass, momentum, and angular momentum at spatial infinity, 
and this excludes all perturbations within the family. Thus, for a family of black holes, the 
thermodynamic stability of the family provides absolutely no information about dynamical 
stability. In particular, this explains why the Schwarzschild black hole can be dynamically 
stable despite the fact that it is thermodynamically unstable: The thermodynamically 
unstable perturbations of Schwarzschild have a nonvanishing linearized ADM mass and 
thus do not "count" for the analysis of dynamical stability. However, the situation is 
very different for a family of black branes. If the corresponding family of black holes is 
thermodynamically unstable — i.e., if we can find a perturbation to another member of 
the black hole family that makes £ < — then we shall show in section 5 that one can 
find a sufficiently long wavelength perturbation of the black brane for which £ < but 
the linearized ADM mass, momentum, and angular momentum vanish. Consequently, the 
thermodynamic instability of the corresponding family of black holes implies dynamical 
instability of the black branes, as conjectured by Gubser and Mitral). 

In this paper, we will restrict attention to black holes (and black branes) in vacuum 
general relativity. However, one of the strengths of our approach is that it is based 
upon the Lagrangian formulation of general relativity and, hence, can be applied to non- 
vacuum general relativity as well as other diffeomorphism covariant theories of gravity. In 

7 In the stable case, the arguments of the previous paragraph strongly suggest that perturbations in 
T must decay to solutions that are perturbations towards stationary black holes. 

8 For the case of spherically symmetric perturbations of static, spherically symmetric spacetimes with 
appropriate matter fields, the linearized constraints can be solved. In that case, our variational principle 
(when generalized to include matter) would reduce to the unconstrained variational principle given in |35] . 

9 Another argument in favor of the Gubser Mitra conjecture based on a "Wick rotation" , applicable 
to certain static black branes, has previously been given by |31j . 
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particular, the notion of canonical energy can be readily generalized to these cases, and it 
appears likely that analogs of properties (l)-(4) above will hold. Furthermore, property 
(5) should hold in non- vacuum general relativity with matter satisfying appropriate energy 
conditions, whereas the satisfaction of property (5) in other theories of gravity should be 
intimately related to the positivity of energy flux at infinity and the satisfaction of the 
second law of black hole mechanics in these theories. Thus, it is likely that most of our 
stability analysis can be applied to much more general cases. However, we shall defer to a 
subsequent publication the analysis of black hole and black brane stability in non- vacuum 
general relativity and other theories of gravity. 

Finally, we should emphasize that, as stated above, for rotating black holes our stabil- 
ity analysis is restricted to perturbations that are axisymmetric in the planes in which 
is nonvanishing. This restriction appears to be essential since, for rotating black holes, 
one should easily be able to find non-axisymmetric perturbations for which £ < and for 
which the flux of canonical energy through J# y+ (or J? + if one instead uses the horizon 
Killing field to define £) is negative. These properties should preclude the possibility of 
making either a stability or an instability argument of the type described above. 

In the next section, we impose our horizon gauge conditions, define canonical en- 
ergy, and establish many of its properties. In section 3 we prove our results on the 
flux of canonical energy at infinity and through the horizon. Our stability /instability 
arguments — including the construction of the Hilbert space V and the domain, T, of 
£ — are given in section 4. In section 5, we present a proof that the thermodynamic in- 
stability of a family of black hole implies the existence of a long wavelength perturbation 
of the corresponding black branes with negative £ that has vanishing linearized ADM 
mass, momentum, and angular momentum. This implies instability of the black brane, 
as conjectured by Gubser and Mitra. In section 6, we shall show that for a family of 
black holes, satisfaction of the local Penrose inequality ([3]) is equivalent to positivity of £ 
on perturbations with vanishing linearized ADM mass, momentum, and angular momen- 
tum, and, consequently, for black holes, the satisfaction of the local Penrose inequality is 
equivalent to linear stability. 

2 Canonical Energy 

2.1 Gauge choice near the horizon 

In this section, we consider one-parameter families of metrics g a b{^) (not necessarily sta- 
tionary) in D > 4 dimensions, that solve the field equations G a b(\) = 0, are jointly 
smooth in the parameter A G M. and spacetime point, and which, for A = 0, reduce to a 
stationary, asymptotically flat black hole spacetime g a b = fi'ab(O), with bifurcate Killing 
horizon Jif, bifurcation surface B C J^ 5 , and surface gravity k > 0. We also will consider 
corresponding families of black brane spacetimes, but the extension of our considerations 
to black brane spacetimes is straightforward, so we will restrict attention to black holes 
until section 5. 

We shall require g a b{X) to be asymptotically flat in the sense that there exist asymp- 
totically Minkowskian coordinates x M such that (^(A) differs from rj^ by 0(p~( D ~ 3 )) as 
p — > oo — where p = (x\ H — • + x 2 D _ 1 )^ — and iVth derivatives of g^X) are 0(p~ ( - D ~ 3+N ' ) ). 
In D = 4, we also require the Regge-Teitelboim parity conditions (see appendix A of [32J) 
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to hold. In section 3, we will, in addition, impose asymptotic flatness conditions on the 
metric at 

We shall now impose certain gauge conditions on the metric near the horizon that will 
be used in the remainder of this paper. These conditions involve no loss of generality. Let 
^ denote the black hole event horizon for the metric g a b{^ — 0), and = t ^'n/ ± (J^ =F ) 
the future and past component. When A > 0, the surface lying at the same coordinate 
location as Jrff + need not coincide with the actual event horizon of g a b{^) — indeed, there 
may not even be an event horizon when A > 0. It will be very useful to require that our 
coordinates be chosen so that, when there is a horizon for A > 0, the surface lying at 
the same coordinate location as J^f + (we could equally choose J4?~) continues to coincide 
with the actual event horizon, at least to first order in A. Since the event horizon is a 
global concept and cannot easily be located exactly in any spacetime with A > 0, this 
might appear to be a hopelessly difficult task. However, as we shall see, the event horizon 
can be easily located to first order in A, and the desired gauge condition can be imposed 
as follows: 

First, by using Gaussian null coordinates, we can ensure that for all A, the surface 
located at the same coordinates as Ji? + is a null surface. Specifically, we may assume 
without loss of generality, that, near Jtf y+ , each metric in our family can be written in the 
form 



where r, u are functions defined in a neighborhood of that are independent of A. The 
location of for any value of A is given by r = 0, and on M ?+ , u is an affine parameter 
for any metric in our familyf^. The tensor fields (3 a , [i a b are orthogonal, for all A, to the 
normal bundle of the (D — 2)-dimensional surfaces B(u,r), which are the joint level sets 
of u, r. That is, if 



then we have /i a ;,n a = [i a bl a = (3 a n a = f3 a l a = 0. For all A, the integral curves of l a are 
null geodesies with affine parameter r. For more details on Gaussian null coordinates, see 
e.g. [21] • We assume that u is chosen such that B = B(0,0) is the bifurcation surface 
of the original metric <? o &(0). The location of the past horizon at A = is then given by 
u = 0. The above form of the metric near ^ + is merely a choice of gauge, which we will 
impose for the rest of this paper. It does not use any information coming from Einstein's 
equation, nor any symmetry assumptions. 

Although the above gauge choice requires J4? + to be a null surface for all A, it does 
not impose the desired condition that J4? + coincide with the actual (future) event horizon 
to first order in A. To impose this condition, we note that the future event horizon is 
uniquely characterized as being the outermost outgoing null surface with compact cross- 
sections whose expansion $ vanishes asymptotically at late times. However, the linearized 
Raychaudhuri equation implies that, to first order in A, the expansion must be constant 
along the horizon. Thus, the perturbed expansion 5$ must vanish on the horizon at all 
times, and, thus, must vanish on any cross-section. Since the requirement that J% y+ is 
null has already be imposed, we can ensure that J% y+ coincides with the event horizon to 

10 Note that in [21], u was chosen to be a Killing parameter rather than an affine parameter, resulting 
in some differences in the metric form in the case of a non-degenerate horizon. 



g ab (X) = 2 V (a w [V b) r - r 2 a(A) V fe) w - r A,) (A)] + // a6 (A) 



(5) 




(6) 
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first order in A by simply imposing the requirement that its perturbed expansion vanishes 
on any one cross-section, such as B. Thus, our desired additional gauge condition is [36] 



8ti\ B = . 



(7) 



We now wish to argue that this condition can always be imposed. Let (p s be the 1- 
parameter family of diffeomorphisms generated by the vector field X a = (l/v2) fl a , 
where / is some smooth function. Let B s = 4> S [B] be the deformed surface, and let 
be the expansion [of the unperturbed metric g a b = g a b(0)] of B s in the outward going 
future null directions k a , normalized by k a l a = 1. The vacuum Einstein equations allow 
one to show that the following equation holds on B [TJ [15] : 



where R(fi) is the curvature scalar of and D a is the derivative operator of /z Q &. Since 
the background spacetime is stationary, we have ^ n fi a b = on Jf? + . As discussed in [TJ, 
although C is clearly not self-adjoint, it admits a strictly positive principal eigenfunction 
■ip. As argued in [TJ [15], we must have Ai > — otherwise there would be trapped surfaces 
in the exterior, which is not possible. We now show that Ai is strictly positive, Ai > 0. 
(For another argument leading to the same conclusion in the more general context of "non- 
evolving horizons", see Prop. 3 of [25] •) To do so, we use the fact that in the background 
metric, B is the bifurcation surface of a bifurcate Killing horizon with horizon Killing 
field K a . We claim first that in our Gaussian null coordinates, K a takes the remarkably 
simple form 



where n is the surface gravity of the black hole. To prove this, we note that the standard 
relation between the affine parameter u on Jif + and the flow parameter of the Killing 
field K a on Ji? + implies that on Ji? + (i.e., when r = 0) we have K a = nud/du. It 
follows that on J$? + , flow under K a by parameter a maps u to uexp(Ka) and n a to 
n a exp(fta). Since l a on is uniquely determined up to scale by orthogonality to (du) a , 
it follows that flow under K a maps l a to Z a exp(— Ka) on '. Since flow under K a maps 
null geodesies to null geodesies, it follows that flow under K a by parameter a acts the 
Gaussian null coordinates as follows throughout the domain of validity of the coordinates: 
u — > uexp(na), r — > r exp(— na), and x A — > x A , where x A denotes the remaining Gaussian 
null coordinates. The infinitesimal form of this transformation is just (Q, as we desired 
to show. 

Now consider the above 1-parameter family of surfaces B s with / = if). For s±, S2 > 
0, the surface B S2 is obtained from the surface B Sl by Killing transport with Killing 
parameter -log(s2/si)- Furthermore, k a = (1/r) k a is Lie-transported by K a . It follows 
immediately that the expansion d(s) with respect to k a is independent of s. From this, 
it follows that the expansion i?(s) with respect to k a is given by i?(s) = s$(l), where $(1) 
denotes the expansion of the surface B\. But we have T"i?(0) = Consequently, if 

Ai = 0, we would have $(1) = 0, i.e. B\ would be a MOTS lying outside the black hole, 
which is impossible. This establishes that Ai > 0. 



d 
ds 



■#(0) = -D a DJ + f3 a D a f + 1 R{») - l - (3 a f3 a + D a f3 a - l - {^ ab )^ ab f 
■ C(f), " ' . 



8 




(9) 
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Thus, we have shown that the stability operator C on B, defined by (jHJ), has a strictly 
positive principal eigenvalue. But the principal eigenvalue of the adjoint C* is equal to 
that of C (see lemma 4.1 of [I]), so C* cannot have a kernel. This implies that we can 
uniquely solve the equation C(f) = j for any smooth j. Solving this equation with 
j = — 5$\b, we find that the 1-parameter family of metrics <p* x g a b{\) will have d\ B = 
to (zeroth and-) first order in A. Thus, we have proven that for any metric perturbation 
7 a b = -j^g a b(0), we can apply a first order gauge transformation 7^ — > 7^ + S£x9ab of the 
form X a = (1/V2) fl a that imposes eq. ([7]). 

We are still free, without destroying our previous gauge choices, to apply a 1-parameter 
family of diffeomorphisms (j) s leaving r, u invariant. Such a diffeomorphism is locally (near 
uniquely determined by its action on the bifurcation surface B, and at the linearized 
level corresponds to a gauge vector field £ a tangent to B satisfying Jzf^r = = S£(U. It is 
possible to use this freedom to impose the additional gauge condition 

H ab 5[i a b\B = const. , (10) 

or, equivalently, 

ft A 

Se\ B = — e\ B (11) 

where e denotes the background volume form on B and S A denotes the perturbed area of 
B. To prove this, we note that Sfj, a b changes under a linear gauge transformation of the 
kind we just described by 2D( a £ a ), where D a denotes the derivative operator on B. Thus, 
to impose the desired gauge, we must solve 2D a ^ a = —fi ab Sfi a b + const. The necessary 
and sufficient condition for obtaining a solution to this equation is that the right side 
integrates to zero over B, which can be achieved by an appropriate choice of the constant. 
It is easily seen that f lTUj) is equivalent to f lTT]) . 

We summarize our gauge choices in the following lemma: 

Lemma 1. Let (^#, g a b) a stationary black hole spacetime satisfying the vacuum Einstein 
equations, and let <? a &(A) a 1-parameter family of solutions perturbing g ab = g ab (0). Then 
near C the metric can be brought into the form (jS]) with 5i9\b = (and, hence, 
the perturbed expansion vanishes on all of Jif + ). Furthermore, the gauge condition fill I) 
can be imposed at B. 

Remark: The gauge freedom that remains on a perturbation 7^ after the imposition of 
the gauge conditions of lemma [TJ is 7^ — > ^ a b + J^?x9ab where X a is a smooth vector field 
that is an asymptotic symmetry near infinity, is tangent to satisfies J^ n X a \^+ = fn a 
with n a \7 a f = 0, and satisfies (n ab \7 a X b )\ B = 0. 

Finally, the black hole background metric g a &(0) is stationary, with Killing field t a . By 
the rigidity theorem, we know that t a can be decomposed as 

N 

t a = K a -^^A a , (12) 

A=l 

where, for A = 0, K a is a Killing field normal to ip\ are Killing fields with closed 
orbits, and Qa are constants called the "angular velocities of the horizon" . On the future 
horizon, we have from (J9j) 

^="»(U- (i3) 
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In our analysis, we will later restrict attention to the axisymmetric case, although we do 
not make any such restriction now. More precisely, we will assume later that for each A 
appearing in the sum in ffT2l . there is a corresponding Killing field of g a b{X) for 

A > 0. Without loss of generality, we may assume that = ip% i.e., it is compatible 

with the gauge choices we have made above to additionally require that our gauge is 
chosen so that the do not vary with A. For A > 0, g a b{X) is not, in general, stationary, 
so there are no corresponding Killing fields t a (A) or K a (X). Nevertheless, it will be useful 
to consider the A-independent vector fields t a and K a in the spacetime (^,g a b{X)), and 
we will do so in the next subsection. However, we caution the reader that, in the event 
that g a &(A) is a stationary black hole, with Killing fields t a (X) and K a (X), we may not 
assume that t a (A) = t a or K a (X) = K a . Indeed, if g a b(A) is a stationary black hole with 
surface gravity or horizon angular velocities differing from that of 5W>(0), it would be 
manifestly inconsistent with eqs. f|T2|) and ffT3]) for t a (A) = t a and K a (X) = K a . Thus, we 
will have to exercise care at the end of section 2.3 when identifying perturbations toward 
stationary black holes, since we cannot assume that these perturbations are manifestly 
stationary in our gauge. 



2.2 First variation formulas 

In this subsection, we will develop the machinery needed to define the canonical energy 
and establish its properties. 

The Lagrangian for vacuum general relativity is 

L ai ...a D — ^g^-R £-a\...a D i (14) 

where R is the Ricci scalar and e ai ,,, aD is the positively oriented^ volume form of a D- 
dimensional spacetime with metric g a b. We will consider the variations of L and other 
quantities under smooth, one-parameter variations of the metric g a ft(A), and, occasionally, 
under two-parameter variations 5W>(Ai, A2). We will write 

Sg a b = jx9ab\\=o , SL = ^£|a=o , S 2 L = £?L\ x =o etc. (15) 

for the derivative (s) at A = 0. In cases where our variational formulas hold for all A, we 
will use the notation ll d/d\ n rather than "5." 
The variation of L can be written as 

±m = E(jg)-±g + Mto±g), (16) 



where E = are the field equationd_f| and 6 corresponds to the boundary term that would 
arise if the variation were performed under an integral sign, namely 

$ai...a d _! — Yq^ V 6ca l— a d-l > (17) 



11 We choose orientations in this paper is as follows: The orientation of ^# is defined by declar- 
ing the Gaussian null coordinates (r, u, x 1 , . . . , x D ~ r ) to be right-handed. The orientation of E is de- 
fined by declaring (r, x 1 ,... ,x D ~ r ) to be left-handed, and the orientation of B is defined by declaring 
(x 1 , . . . , x D_1 ) to be right-handed. 

12 More precisely E ab Cl ... CD = -^e Cl ... CD (R ab - \Rg ab ) 



12 



Uai...a D -i ~ i(i„ w e cai...a D -n (21) 



where 

v a = g ac g bd (V d ±g bc - W c f x , 

The symplectic current (D — l)-form u is defined as 

"Ok ii9, £- 2 g) = atr%; £- 2 g) - £- 2 6(g; &g). (19) 

Thus, oj depends on the unperturbed metric and a pair of perturbations gfj#)j an d 

it is antisymmetric in af~sO- If both perturbations satisfy the linearized equations 

of motion, then it follows by taking a second, antisymmetrized variation of (fl6|) that u is 
closed, 

dco = . (20) 

Concretely, 71,72) is given by 

1 

16-7T 

where 

w a = pabcdef^ ^ ^ ^ _ ^ ^ ^ ^ ^ 

and 

pabcdef = g ae g fb g cd _ \ g ^ g be g fc _ l g ^ g cd g ef _ }_ g bc g ae g fd + }_ g bc g ad g ef ^ 

The symplectic form Ws(g] 71, 72) is obtained by integrating u over a Cauchy surface 

^(fi-; 71,72) = J u(g; j h j 2 ) . (24) 

In the case mostly considered in this paper — where g ab = g a b(0) is a stationary black 
hole with bifurcate Killing horizon (see the previous subsection) — we take S to be a 
Cauchy surface for the exterior region, i.e., a spacelike slice that goes from the black hole 
bifurcation surface B to spatial infinity. For D > 4 this integral manifestly converges for 
our asymptotic conditions on the metric near spatial infinity (see the second paragraph 
of the previous subsection). For D = 4, the Regge-Teitelboim parity conditions assure 
convergence. If (71,72) both satisfy the linearized field equations around a solution g ab 
to Einstein's equations, then then it follows from du = and our asymptotic conditions 
near spatial infinity that 

Ws/(#;7i,72) = W E (#;7i,7 2 ) (25) 

where £' is another Cauchy surface extending from B to spatial infinity. 

We can bring the symplectic form into a more recognizable form by performing a space- 
time split of the spacetime metric into the canonically conjugate variables (h a b,p ab ), where 
h a b is the induced metric on E and 

v ah = hh{K ah -h ah K), (26) 

with K ab the extrinsic curvature of S. Our asymptotic conditions near spatial infinity 
imply t hatful 

h ab = 5 ab + 0(p-^) , p ab = 0{p-V>-*>) , (27) 



13 Here it is assumed that £ approaches a surface of constant a; in our asymptotically Minkowskian 
coordinates. 
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with the N-th spatial derivatives of h ab and p falling off faster by a factor of p . Note 
that we have followed the standard practice of incorporating the volume element of E 
into the definition of p ab , as is necessary in order that (h ab ,p ab ) be canonically conjugate. 
Thus, p ab is a density of weight |. Our integrals below involving p ab will be taken with 
respect to a fixed, nondynamical (e.g., coordinate) volume element e$ on E, such that the 
volume element v a e abl ... bD1 (where v a is the future directed unit normal to E) associated 
with hab, is equal to hh (e )b 1 ...b D _ 1 - A perturbation 5g ab satisfying the linearized Einstein 
equations corresponds to a pair (Sh ab , 5p ab ) satisfying the linearized constraints and the 
linearized Hamiltonian equations of motion. In terms of (Sh ab , 5p ab ) the symplectic form 
as defined by eq. (124"|) is simply0 

Wx(g; S l9 , 6 2 g) = [ (5 1 h ab 5 2 p ab - 6 2 h ab 6 lP ab ) . (28) 

167T J s 

Here, and in the following, we omit the reference volume form e$. 

For an arbitrary vector field X a , the associated Noether current J x is defined by 

J x (g) = e{g;& x g)-i x L(g) (29) 

where i x denotes contraction of X a into the first index of a differential form such as 
L ai ...a D - A simple calculation [21] shows that the first variation of J x satisfies 

±J X = -i x (E ■ ±g) + u(g; &g, Jf x g) + d[i x 0(g, ±g)} , (30) 

where, in this formula, it has not been assumed that g ab satisfies the field equations nor 
that dg a b/d\ satisfies the linearized field equations. Furthermore, it can be shown [25] 
that J x can be written in the form 

Jx = C x + dQ x , (31) 

where C x = C a X a are the constraints of the theory |35j and Q x is the Noether charge. 
In the case of general relativity in vacuum, eq. ffl4|) . we have 

(Cx)oi...or,_i = ^ X a G b e bai ... aD _ 1 (32) 

and 

(Qx) ai ...a D _ 2 = _ Y67T ^ bXcebC ^- a D-2 ■ ( 33 ) 

Combining eqs. fl3"0l and fl3Tl) . we obtain 

u(g; &g, JZ x g) = i x {E{g) ■ ±g) + ±C x {g) + d [±Q x {g) - i x 6(g; ±g)] (34) 

It should be emphasized that eq. (134"|) is an identity that holds for arbitrary X a and g ab (X). 

We now impose the condition that our one-parameter family g a b{X) is composed of 
asymptotically flat solutions to Einstein's equation, so E = C a = 0. We also require that 

14 The pullback of (pH) to a spacelike hypersurface was computed in eq. (4.14) of [2]. There is an 
additional "boundary term" present in their expression for , which is easily seen to vanish on account 
of our boundary conditions at B and our asymptotic conditions at infinity. 
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X a be an asymptotic symmetry near infinity, i.e., we require ££x9ab to satisfy the same 
asymptotic conditions at infinity as 5g a b, so that X a suitably approaches a Killing field of 
the Minkowskian background as p — > oo. Integrating eq. fl34l) over a compact region, K, 
of E, extending from B to a surface, S, in the asymptotic region, we obtain, for all A 

f u;(g; J? x g) = [ [±Qx{g) - &0] - / [&Qx(g) ~ i X 9(g; &g)) . (35) 

Jk Js Jb 

Taking the limit as S — > oo, we obtain 

W^g- ±g, J? x g) = [ [f x Q x {g) - &g)} - [ i£Qx(g) - i x 0(g; &g)) , (36) 

Joo Jb 

where fJ±Qx(jj) - ix9(g; ±g)] = f s [£;Qx(g) ~ p)]. For £> > 4 or for 

the case where X is an asymptotic translation in D = 4, this limit exists independently 
of how S approaches infinity because the integral defining Ws{g; 4tg, ££xg) converges 
absolutely by virtue of our assumed asymptotic conditions. For asymptotic boosts and 
rotations in D = 4, the limit must be taken where S approaches a two-sphere asymptoti- 
cally. 

The formula 

£ H x= [ [&Qx(g)-ixO(<r,fy)], (37) 

J oo 

can be shown to define [39] a conserved quantity H X - If X a — > (d/dt) a , then Hx is 
equal to the ADM- mass M; if X a — > (d/dx l ) a (asymptotic translation), then — ffx is 
equaf^l to the linear ADM-momentum Pf, if X a — > x l {d/dx^) a — x J (d / dx l ) a (asymptotic 
rotation), then — Hx is equal to the ADM-angular momentum Jr^-i in the n-plane^l; if 
X a — > t(d/dx l ) a + x l (d/dt) a (asymptotic boost), then H x is equal to the ADM-center of 
mass, C{. Together, the \D(D + 1) quantities (M, P, J, C) will be referred to as 'ADM 
conserved quantities" associated with the spacetime (^,g a b)- 

Now consider the case X a = K a , where K a is the horizon Killing field of 5w,(0). Then 
K a takes the coordinate form ffT31) at the horizon for all A. In particular, K a \s = and, 
thus, ikO = on B. On the other hand, using ([5]), we obtain 

/ QMX)) = £-A{\) (38) 
Jb 8tt 

where A(\) is the area of B in the metric g a b{X). Thus, we obtain 

S&K{g)-iK9{g ] f x g)\ = ^, (39) 
which holds for the metric g a b = g a bW- By eqs. (IT2|) and ( 157)) . we have 

Jj&QM - i K e(g; ±g)] = ±M - £ n A ±J A . (40) 

15 The map X — > Hx is a linear functional and thus is a "covector" . To get the ZJ-dimensional energy- 
momentum vector, we have to "raise the index" with the Minkowskian metric. This accounts for the 
relative minus sign between M as compared with Pi and Ju^ in their definitions in terms of Hx ■ 

16 The notation J a used above refers to N distinguished mutually orthogonal 2-planes defined by the 
N rotational Killing fields i/j a . 
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Combining eqs. fl36|) . fl39|) . and f )40|) . we obtain 



W^± g ,Jf K9 )^ x M-^n A p A -±±A. (4!) 

This equation holds for all A, with k and Q A fixed constants, equal, respectively, to the 
surface gravity and angular velocities of the stationary black hole metric «7 a &(0). Since K a 
is a Killing field of g a &(0), when A = we have 5£kQ = and eq. (JIB reduces to the first 
law of black hole mechanics 

= 5M - V s n A 5J A - —SA , (42) 

A 

where we remind the reader that our use of the notation "5" here (rather than "d/dX") 
indicates that this equation holds only at A = 0. 

Finally, we consider the gauge dependence of W^. It is clear that W^(g;Sig,S2g) 
will be gauge invariant if and only if W^(g; Sg,^g) = for all allowed Sg a b and £ a . It 
follows immediately from integrating f !34p over £ (with E — 0) that if 5q n .h satisfies the 
linearized constraints, then W^(g; Sg, =£%?) = for all £ a of compact supporO away from 
B. Thus, for solutions, W^(g; S±g, S2g) is gauge invariant for gauge transformations of 
compact support. However, if £ a is a smooth vector field such that ^g a b satisfies our 
asymptotic conditions at infinity (i.e., £ a is an asymptotic symmetry at infinity) and our 
boundary conditions at B, then W^(g;5g,^g) need not vanish for all Sg that satisfy 
our asymptotic conditions at infinity and our boundary conditions at B. The following 
lemma characterizes the extent to which Wy, is gauge invariant. This lemma will play a 
critical role in our analysis. 

Lemma 2. Let 5g a b be a solution to the linearized Einstein equations around our sta- 
tionary black hole background g a b satisfying our asymptotic flatness conditions and our 
gauge conditions ([7]) and ffTTj) at B. Suppose in addition that 5A = (so that, by 
( TTTj) . we have 5e\b = 0) and that 5Hx = for some asymptotic symmetry X a . Then 
W / s(5 1 ; Sg, <=£%?) = for all smooth £ a such that (i) £ a |# is tangent to the generators of 
Jrf? and (ii) £ a approaches a multiple of X a as p — > oo. Conversely, if 5g a b is smooth and 
asymptotically flat and if W-s(g; Sg, Jz^g) = for all such £ a , then Sg a b is a solution to the 
linearized Einstein equation with Si9\b = Se\s = at B and with SHx = 0. 

A proof of this lemma is given in appendix A. 

2.3 Second variations and canonical energy 

Returning to eq. ( |4T|) . we write out the A-dependence explicitly as 

v 7 A 



17 Conversely, if Wz(g; Sg,Jf^g) — for all £ a of compact support away from B, then 5g a b satisfies the 
linearized constraints. These facts, which are manifestations of the familiar statement that in general 
relativity "the constraints generate gauge transformations," will be exploited heavily in subsection 4.1 
below. 
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where it should be emphasized that k and VLa do not depend on A. We now take a A- 
derivative of this equation, and then set A = 0. Using the fact that Jz?Kg a &(0) = 0, we 
get 



2 



Wv{g;i,^Ki) = -^M{\) 



2 



A=0 A aA 



2 



k d 

A=0 ^71 ^ 2 



A(X) 



(44) 



A=0 



where 7 & = dg a b/d\\ x=0 = Sg a b(0) is the first order perturbation of the black hole space- 
time. Writing K a = t a + Y2a ipA a , the left side may be written as 

Wx (g; 7 , S?m) = £ + ^ Oi We {g; 7, ^ A l) , (45) 

A 

where we have defined the canonical energy by 

£ = W E G7;7,J2f t7 ) • (46) 

Since t a is a Killing field of the background, S£ t commutes with the linearized Einstein 
operator, so J^tlab satisfies the linearized field equations whenever j a b does, ^ftjab also 
satisfies our asymptotic conditions at infinity (with faster fall-off) and our boundary 
conditions at B. Hence, it follows from the conservation of symplectic product that £ is 
conserved for all solutions 7^ in the sense that it takes the same value if evaluated on 
another Cauchy surface £' extending from B to spatial infinity. We may rewrite (1441) as 



£ = - ^A W s (g; 7, ^7) + 5 2 M - ^ Q A 5 2 J A - ^6 2 A , (47) 



where we have written 5 2 M = d 2 M/d\ 2 \\=o, etc. Note that £ and Wt, (g; 7, =^/>a7) 
depend only upon 7^ = 5g a b, whereas 6 2 M, 8 2 Ja, and 5 2 A depend upon the second order 
perturbation 5 2 g a b = d 2 g a b/d\ 2 \\ =0 . However, the difference between two second order 
perturbations will be a homogeneous solution of the first order perturbation equations, 
so its net contribution to (14T1) will vanish by the first law of black hole mechanics (1421) . 
We summarize our main result thus far as follows: 

Proposition 1. (Second variation formula): Let g a &(A)) be a one-parameter family 
of smooth, asymptotically flat vacuum solutions such that g a b{0) is a stationary black hole 
spacetime with bifurcate Killing horizon. Suppose that the gauge condition ([5]) has been 
imposed near the horizon J4? + . Then the canonical energy fT46|) of 'jab satisfies (14T|) . where 
Qa and k are the angular velocities and surface gravity of the horizon of the background 
solution g a &(0). 

It is useful to view the canonical energy as a bilinear form on perturbations, defined 

by 

£(71,72) = W E (g; 71, ^72) , (48) 

so that the canonical energy originally defined by (146]) is given by £ = £(7,7). Although 
it is not manifest from its definition, it is easily seen that £(71, 72) is symmetric in (71, 72): 

Proposition 2. Let 71 and 72 be smooth, asymptotically flat solutions to the linearized 
field equations on the stationary black hole background g a b{0)- Then, we have 

£(71,72) = £(72,71) • (49) 
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Proof: Since the background metric g is stationary, it follows immediately from the anti- 
symmetry of the symplectic current u that 



71, JSftTfc) - 72, -2t7i) = -S^fo; 71, 72) • (50) 

Now we use the standard identity S£x = ix d + d %x for the Lie-derivative acting on forms, 
together with the fact that do; = when 71 and 72 are linearized solutions. Integrating 
over S, we obtain 

£(71,72) -£(72,71) = / *M#;7i,72)- / i t uj{g; ji, 72) ■ (51) 

Joo J B 

However, the contribution from infinity vanishes on account of our asymptotic condi- 
tions, whereas the contribution from B vanishes because t a is tangent to B. Thus, £ is 
symmetric. □ 

In our analysis of the next sections, we will be concerned with the positivity properties 
of £ as well as fluxes of £ through the horizon and at null infinity. Unfortunately, for a 
rotating black hole, the term YIiA^aWy, (g', 7, -^ a t) wm s P°il an y positivity properties 
of £ . Furthermore, the fact that, for the case of a rotating black hole, t a is spacelike at the 
horizon willpreclude the possibility of having a positive flux of canonical energy through 
the horizon! 18 !. Fortunately, both difficulties can be resolved by restricting consideration 
to axisymmetric perturbations 

jsV? = o, (52) 

and we shall impose this restriction later in our analysis. For axisymmetric perturbations, 
( H7|) reduces to 

S = 5 2 M-Y n A 5 2 J A - — 5 2 A. (53) 

A 

An important property of axisymmetric perturbations is given by the following propo- 
sition: 

Proposition 3. The canonical energy is gauge invariant with respect to gauge trans- 
formations that preserve the gauge conditions of subsection 2.1 when restricted to the 
space of axisymmetric perturbations that satisfy the linearized field equations and have 
vanishing linearized ADM linear momenta, 5Pi = 0, and vanishing change of area, 5 A = 0. 

Proof: The gauge transformations that preserve the gauge conditions of subsection 2.1 
and axisymmetry are 7^ — > 7^ + J£x9ab where X a satisfies the conditions in the remark 
below lemma[T]of subsection 2.1 and, in addition, satisfies ££^ A X a = everywhere. Since, 
by Proposition 2, £(71,72) is symmetric in (71,72), gauge invariance will be proven if we 
can show that 

£(%S?x9) = W s (j,^ x g) = (54) 



18 We could have defined a canonical energy Ek with respect to the horizon Killing field K a instead of 
t a , in which case the term O^Ws (g; 7, Jzf^ A 7) would not have appeared in the expression for £ k and 
the flux through the horizon would be positive. However, in the absence of axisymmetry, the spacelike 
character of K a near infinity for a rotating black hole will similarly spoil the positivity properties of £k 
and preclude the possibility of its having a positive flux through null infinity 
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for all solutions 7 a & for which 5Pi = and 5A = 0. However, we have 

^t^XQab — \J£t^X — ^X^t)9ab = ^[t,X]9ab , (55) 

since S£ t g a b — 0. Since X a is an asymptotic symmetry, [t, X] a is an asymptotic spatial 
translation. Furthermore, since t a = K a — Y^^a^a an d [^A, X] a = S£^ A X a = 0, it 
follows that [t, X] a \B is normal to Jif + at B. The proposition now follows immediately 
from lemma [2] applied to £ a = [t, X] a . □ 

In our later arguments below we will need to consider perturbations towards other 
stationary black holes. We now explain the subtleties involved in this notion and give a 
precise definition of it. It might appear to be obvious from its definition (j4"6"|) that £ must 
vanish for stationary perturbations. However, this is not the case because, as already 
indicated at the end of subsection 2.1, we cannot assume that the stationary Killing field 
t a (X) of the varied spacetime coincides with t a . Indeed, the example of the one-parameter 
family of Schwarzschild metrics of mass M = Mq + aX in 4-dimensions illustrates that 
the canonical energy need not vanish for perturbations to other stationary black holes. 
For this family, we have S 2 M = 0, but A(X) = 16tt(M + aA) 2 , so 5 2 A = 32vra 2 . Thus, 
by dSSD, we have (using k = 1/4M ) 

£ = -a 2 /M < (56) 

To see why this result is compatible with (1461) and to find the appropriate conditions 
to express the notion that a perturbation 7^ represents a perturbation towards another 
stationary black hole, let g a b{X) be a one-parameter family of stationary black holes, with 
stationary Killing field t a (X), horizon Killing field K a (X), and axial Killing fields ipA a (X) 
associated with a nonvanishing angular velocity of the horizon. As already stated at the 
end of subsection 2.1, it is compatible with our other gauge conditions to assume that 
the axial Killing fields ipA a {X) are independent of A, so that ^"(A) = ipA a - We may also 
assume that near infinity, we have t a (X) = t a . However, since K a takes the form ffTS]) on 
we must have K a (X) 7^ K a near Jif + if the surface gravity, k(X), of g a b{X) differs 
from k. Nevertheless, we may assume that K,K a (X) = n(X)K a near J4? + . If we make this 
choice, then, near J4? + we have 

t« (A) = <^ K a _ J2 n A (x)^ A a , (57) 

A 

so t a (X) must differ from t a near Jif + if there is any change in the surface gravity or 
angular velocities of the horizon. 

Thus, if g a b{X) is a one parameter family of stationary black holes, then we can choose 
a gauge compatible with our previous choices such that ^"(A) = ipA a , t a (X) = t a near 
infinity, and t a (X) is given by (15?]) near J4? + . The perturbation 7^ = dg a b/dX\\=o must 
therefore satisfy ~^ A 7 a 6 = for all ipA a appearing in (157|) and 

^tlab + ^ab = (58) 

where, near infinity, we have St a = whereas near J4? + , we have 

St a = c t a + J2 b A ^A a (59) 

A 



19 



where c and 6 a are constants. Thus, near 5t a is of the form of a Killing field of the 
background metriJ^l. 

Thus, rather than having to vanish, we see that J^tlab must take the form Jzfyg a f> 
where Y a vanishes near infinity and is a linear combination of Killing fields ( 15 9p near . 
Equivalently, writing Z a = Y a — ct a — J2a ^a^a 1 and noting that J^yQab = ^z9ab-, we 
see that a perturbation towards a stationary black hole can be put in a gauge compatible 
with our gauge conditions such that 

^tlab = ^Z9ab , (60) 

where Z a vanishes near Jif + and near infinity takes the form 

Z a = ct a + J2b A ^A a , (61) 

A 

and where, furthermore, =Sf^ A 7 a fe = for all ipA a appearing in ( TT2"|) and/or ( 16"T|) . In a 
general gauge satisfying our gauge conditions, we have 

^ tlab = ^9ab + &t&X9ab , (62) 

where X a satisfies the conditions of the remark below lemma [TJ Thus, taking account of 
f )55|) and the remarks below that equation, we have have motivated the following definition: 



Definition 2.1. A smooth, asymptotically flat, axisymmetric solution 7^ of the linearized 
field equations is said to be a perturbation towards a stationary black hole if 

££ tlah = Sftfa . (63) 

where £ a |e is normal to Jif + and near infinity takes the form 



A o" \ / 



(64) 



Note that although the linearization of any one-parameter family of stationary black 
holes must satisfy definition 2.1, our definition does not require that there actually exist 
a one-parameter family g a b{^) of stationary black holes corresponding to jab- In any 
case, the key point about perturbations towards a stationary black hole is that they are 
stationary and thus are benign with regard to linear stability of the background black 
hole. An important property of perturbations towards a stationary black holes is the 
following: 

Proposition 4. Let 71 be a perturbation for which 5M = 5Ja = 5Pi = (and hence, 
by the first law of black hole mechanics 5 A = 0), and let 72 be a perturbation towards a 
stationary black hole. Then £(71,72) = 0. 

19 The sum on the right side of f|59[) is allowed to include rotational Killing fields ipA a of the background 
that may not have appeared in (TT21) because ^^(A = 0) =0. 
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Proof: We have 

£(71,72) = W^G?;7i,^7 2 ) = Wv(g;-yi,Sftg) = 0, (65) 

where the second equality merely substitutes the definition of a perturbations towards a 
stationary black hole and the last equality is an immediate consequence of lemma [2] □ 

In subsection 4.1 below, we will considerably strengthen this result by showing that, 
when considered as a quadratic form on a suitable space of smooth, axisymmetric solutions 
to the linearized field equations for which 5M = 5Ja = 5Pi = 0, the canonical energy 
£ will be degenerate precisely on the perturbations towards stationary black holes. This 
will play a key role in our stability arguments. 



2.4 Evaluation of canonical energy 

In the previous section, we obtained a simple formula for £ in terms of second order 
variations of mass, angular momentum, and area. However, in order to use this formula 
to evaluate £ , we must calculate second order perturbation d 2 g a b/d\ 2 \\ = o, even though 
£ really only depends on the first order perturbation 7^ = dg a b/d\\\ = o. It is useful to 
have a formula for £ that expresses it directly in terms of 7 a t, and/or the initial data, 
(Shab, 5p ab ), for '-fat,. Such a formula can be obtained from the original definition 



£(7i, 72) = W^g- 7l , J2? tT2 ) = / fak* J? t 5 2 p ab - 5 lP ab J? t 5 2 h ab ) , (66) 

by substituting the explicit expressions (12TT) and (1221 . However, to put the right side 
in a more useful form, we will follow a different strategy and return to the fundamental 
identity (I34p . using the fact that the variation is being taken about a solution, so E = 0. 
We obtain at A = 

co(g; Sg, J? x g) = 5C x (g) + d [5Q X - i x 9(g; 6g)} . (67) 

The constraints C x = X a C a are given by eq. (1321) above. In terms of the variables 
(h a b,p ab ), the constraints take the fornix 

C« = (Cf t ) := 1 Hi + h l P €r t t^ V ) - , (68) 



CbK) ' 16tt V -2D b (h-2 Pa b 

where v a denotes the future-directed unit normal to E. The first line of this equation 
corresponds to the Hamiltonian constraint and the second line corresponds to the mo- 
mentum constraint. The linearized constraints 5C a may be viewed as the result of acting 
on (Sh a b, 5p ab ) by a linear operator, £, 

6C a = C fe) (69) 



20 Here and below, we omit writing the factor of the non-dynamical coordinate (D — l)-form eo on E 
in the expression for C a [see the discussion below (|27p]. 
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where C is explicitly given by 



16tt£ 



5h a 
5p a 



( hl{D a DJh c c - D a D b 5h ab + R ab (h)5h ab )+ \ 
h-l{-5h c c p ab p ab + 2 Pab 5p ab + 2p ac p b a 5h bc + 
1 p c c p d d Sh a a - -fjzoP a a 5p b b - Thl 5h abP ab Pc c 



D-2' 



(70) 



\ -2h-2D b (h-2Sp ab ) + D a 5h cbP cb - 2D c 5h abP bc J 



Since £ is a differential operator that maps the pair (8h ab , Sp ab ) consisting of a symmetric 
tensor, Sh ab , and a symmetric tensor density, 5p ab , on £ into a pair (N, N a ) consisting of 
a scalar density and dual vector density on S, its adjoint differential operator, £*, maps a 
pair (iV, A^ 1 ) consisting of a scalar and vector field on £ into a pair (5h ab , 5p ab ) consisting 
of a symmetric tensor density and symmetric tensor on S. In other words, C, C* are maps 



£ : C°°(E, (T*S) W ) © C°°(E, (TE) W © A*) -»■ C 00 ^, A*) © C°°(S, T*£ © A*; 



£* : E) © TE) ->■ (TS) V2 © As) © (T*£ 



(71) 



where A^ is the line bundle of densities of weight |, and V is the symmetric tensor product. 
C* is uniquely determined as a differential operator by the requirement that for all smooth 
(Sh ab , Sp ab ) of compact support on £ \ B and all smooth (N, N a ), we have 



N 
N a 



5h ab 
5p ab 



N 
N a 



C 



Sh ab 

5p ab 



(72) 



Here, the angle brackets denote the natural L 2 -type inner product on the appropriate 
tensor/tensor density pairs, e.g. 



N 
N a 



N 

Nn 



(NN + N a N a ) 



(73) 



One can straightforwardly calculate that C* is given by 



16tt£* 



N 
N a 



( hh (—D c D c Nh ab + D a D b N + R ab (h)N)+ \ 
h-\(-h ab p cd p cd N + 2p { - a c p b ^N + -^h ab p c c p d d N 



p ab p c c N - p ab D c N c + 2 J D c My )c ) 



V 



D-2 

h- X 2{2p ab N 



(74) 



D-2 



h ab p c c N) + 2D (a N b) J 



Comparing ( 172"]) with the integral of f )67j) over £ with (Sh, dp) of compact support on 
£ \ B, we see that 



8p 



Here by C*(X), we mean C*(N,N a ), where 

X a = Nu a + N a 



(75) 



(76) 
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is the usual decomposition of a vector into its lapse and shift. Since neither side of 
contains derivatives of (5h,Sp), we can solve for JZ'xd, obtaining 



J?: 



x 



16tt aC*{X) , a 





_ h ac h bd 



h ac hbd 




(77) 



which corresponds to the usual ADM evolution equations. Note that — in view of (177)1 — 
(175]) continues to hold even when (Sh, Sp) is not of compact support on £ \ B, provided 
only that the integrals defining the left and right sides converge. Note also that when 
(5h, 5p) is not of compact support on £ \ B, integration of ( 1dT|) over E yields at A = 



C*(X) 



5h 
Sp 



X 



5h 
dp 



+ 



[5Q x (g)-ixe(g;5g)\ 
[6Qx(g)-ixe(g;Sg)} 



(7? 



B 



To derive the desired formula for £(71, 72), we now consider a two-parameter family of 
metrics g ab (Xx, A 2 ) = g ab + A^ a6 + A 2 72 ab, where is the background metric and 71 ab 
and 72 a b solve the linearized equations. Although g a b(Xi, X 2 ) is not a solution, it does 
satisfy the field equations to first order in both Ai and A2. We find that the corrections 
to (T75|) away from Ai = A2 = are of the form 



W^[g(X 1 ,X 2 ); 



± i g(0,X 2 ),^ x g(X 1 ,X 2 ) 



C*(X) 



+ 0(X 1 )+0(Xl). (79) 



Taking the derivative of this equation with respect to A 2 , setting Ai = A 2 
also setting X a = t a , we obtain 



£(71, 72) = We (g; 71, ^72) 



0, and then 



JO) 



Ai=A 2 =0 



The right side of this equation can be computed by taking the variation of the right side 
of ()74p . although we shall not explicitly write out the resulting formula here. Note that 
although 5t a = 0, it is not compatible with the gauge conditions we have already imposed 
near the horizon to require 5v a = 0, so variations of the lapse and shift must be taken 
into account in the calculation of (t)\\ 1= \ 2 =o- 

An alternative formula for £ can be derived by noting that, to the desired accuracy, 
one may substitute ( 17H|) on the right side of (179|) . Taking the derivative of the resulting 
equation with respect to A2, and again setting Ai = A2 = as well as X = t, we obtain 



£(71,72) 



d\2 



Ai=A 2 =0 



+ 



B 



(81) 



It is easily seen that the boundary contribution from infinity vanishes. Furthermore, the 
pullback to B of i t -£^9(g; -g^g) vanishes for Ai = A2 = 0, since t a is tangent to B. Thus, 
we obtain 



£(71,72) 



_d_ 

d\ 2 



C 



d 2 
8X18X2 



Q 



t , 



2) 
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where evaluation at Ai = A2 = is understood again. On the other hand, we have 



2 



Ai=A 2 =0 



d\id\ 2 i E 



t a C a (g + X 1 5 1 g + X 2 6 2 g) 



(83) 



Ai=A 2 =0 



Since t a C a = -±G ab v a t b the right side may be computed by taking second order 



variations of the Einstein tensor. This gives 



a=o JB 



dX 



;Qt(g + A7) 



A=0 



Explicit evaluation of the two terms on the right side (see [20J for the first term, and 
use ([2]), (fT2"l) . (I3"3"j) . dSJ) for the second term) gives: 

e = ^ J^(^ cd V a V blcd -2 1 cd V c V [alb)d + ^(V al ^ 

+V d { 1 dc V clab ) - i(V c 7 )V c7 a 6 - 2(V d7 cd - ^V c 7 )V (a7 6) c - ^(trace)) v a t b , 
/x5 <5/x a6 <^° 6 , (85) 



167T 

where 5fi ab is, as before, equal to the pullback of 7^ to B, where "trace" denotes the trace 
of the preceding terms, and where use has been made of the gauge conditions at B in the 
computation of the boundary term. For later use, we also quote the lengthy expression 
for £ in terms of the variables (Sh ab , 5p ab ) and the lapse and shift (N, N a ) of t a , which 
can be obtained by taking a second variation of fl68l) : 



£ 



^ J N^!^R ab (h)Sh c c Sh ab -2R ac (h)Sh ab Sh 



16vr _ 

l - 6h ac D b D b 6h ac + 5h ac D b D a 5h cb - ^ D a (5h bc D a 5h bc ) - ^ D a (6h ab D b 6h c c ) + 
l - D a (6h/D a 6h c c ) + 2 D a (5h a c D b 5h cb ) + D a (5h b c D b 5h ac ) - i D a (5h c c D b Sh ab ) 



1 



h~* |2 5 Pab 5p ab + - Pab p ab (5h a a y - Pab 5p ab 5h c c - 3 p a b p bc 5h/5h ac - 

2 (Spa a ) 2 + 7^ Vc c 5p b b 5h a a + Pd d p ab 6h c c 5h ab + 8 P c b 5h ac 6p ab + 



D - 2 v " 7 D — 2 

p c dP c "^a6^ ab + 2 p ab p dc 5h ac 5h 



D-2 
1 



6d 



D-2 



(p c c Y5h ab 6h 



ab 



1 



2(L> - 2) 



(P 6 b ) 2 (^a a ) 2 



1 Pc c 5p ab 5h ab - — ^ (p ab Sh ab ) 2 - 7 ^ Pab Sp c c 6h ab 



D-2 
1 



167T 



D — 2 ' D-2 

N a [ -2 5p bc D a Sh bc + 4 5p cb D b 5h ac + 2 5h ac D b Sp cb - 

K 



2 p cb 6h ad D b 6h c d + p cb 5h ad D d 5h cb 



16n 



B 



(i* ( Sfi ab 5fi ab - - 6n a a 5n b b 



(86) 
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3 Flux of canonical energy at infinity and the horizon 



As shown in the previous section, canonical energy £(7) is conserved for all solutions 7^ 
of the linearized Einstein equation in the sense that it is independent of the choice of 
Cauchy surface S extending from spatial infinity to B. However, if we evaluate £ on a 
slice S^(t) that extends from a cross-section, C(t), of future null infinity J^ + , (rather than 
from spatial infinity) to a cross section, B(t), of the future horizon, J# y+ , (rather than to 
B) 

£{ 1 ,y{t)) = W y(t) { 1 ^ tl )= [ w(<7; 7 ,Jgf tT ), (87) 

Jy{t) 

then, of course, in general we will find that £(7, ^(t)) 7^ £(7)- In this section we study 
the time evolution of £(7, ^(t)) as we march the slice S^(t) forward in t. We will show 
that, up to boundary terms, £(7, J5^(t)) decreases with t. 

In order to be able to make use of the machinery of null infinity, we will restrict 
consideration in this section to even dimensional spacetimes — J? does not exist for odd 
dimensional spacetimes |23j. However, we do not believe that this is an essential restric- 
tion, i.e., we believe that our results hold in odd dimensions, with a suitable notion of 
asymptotic flatness for that case, see e.g. [37] . 

We work in the "Bondi gauge" for the background metric g abl so the unphysical back- 
ground metric g ab = 0?g ab near takes the form 

g ab = 2\/ (a nV b) u + fx ab + 0(n) (88) 

with u a future directed affine parameter on the null geodesic generators of On 
we have (d/du) a = g ab, V b Q = n a , and jl a b is the unit round metric on the S D ~ 2 cross- 
sections of J^ + , with jl ab n a = = jl ab (d / dVt) a . The asymptotically timelike Killing field 
t a can be extended continuously to a vector field t a on which is proportional to n a , 
where we consequently have 

t a = (t c W c u) n a , (89) 

with t c V c u > and constant on J? . If the initial data for 7^ on S are of compact 
support^, then it has been shown in [52J that there exists a gauged near future null 
infinity in which 7^ is asymptotically flat at future null-infinity in the sense that 

% b := Q-^ 6)/2 lab (90) 

is smooth on J^" 1 ", and that iabg ab = 0(f2),7a & n a = 0(0,). 

It is also assumed throughout this section that 7^ is axisymmetric with respect to all 
of the axial Killing fields appearing in (|12|) 

^7ab = for all A = 1, . . . , iV . (91) 

This restriction is essential because although we will obtain positivity of flux results at 
for canonical energy, we will obtain positivity of flux results at for the analogous 

21 Below, we will consider perturbations whose initial data are in a space T . Elements in this space 
may not have compact support, but it is shown in lemma[5]of Appendix C that the initial data of compact 
support are dense in this space. 

22 In D = 4, this is the Geroch-Xanthopoulos gauge, whereas in even D > 4, it is the transverse-traceless 
gauge, see [22] for further details and other subtle differences between D = 4 and the rest. 
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quantity denned by replacing t a by the horizon Killing field K a . It is only in the presence 
of axisymmetry fl9T|) that these quantities are equal. 

Let 0^12 denote the portion of future null-infinity bounded by cross-sections C(t\) and 
C(t 2 ), and let M[ 2 the portion of the horizon bounded by the cross-sections B{ti) and 
B(t 2 ), see the following figure. 



BH = Jt\ J~{^ + ) 




Clearly 

y(t x ) u jt 12 u j X2 u y{t 2 ) (92) 

bounds a compact region of the conformally compactified spacetime. By conservation of 
canonical energy, we have 

£(7,^2)) = £(7,^(*i)) - w^^a) - w.^^a) . (93) 

Thus, £(7, y{t 2 )) differs from £(7, S*(ti)) by the "flux terms" W^ 12 (j 7 &a) and Wjr 12 (7, ^7), 
corresponding, respectively, to the canonical energy radiated to null infinity and into the 
black hole. 

We now evaluate W^ 12 (7, ^7). Using a formula from [22J based on the Einstein 
equations for the background and perturbations, we have 

WV 12 (7i, 72) = I (7i cdS 2 N cd - 72 cd 5iN cd ) , (94) 

where here and in the following, the natural integration measure on J? + coming from g a ^ 
is understood!! and where indices on tilde tensor fields are always raised and lowered 
with g a b. The Bondi news tensor is defined b}@ 

N cd = ffcitd n-^ 



23 In local coordinates fi, u, xi, . . . ,xd-2 adapted to (j88|) . this would be p,^dudxi ■ ■ ■ dxp-2- 

24 The form of the second term depends in general on the the conformal factor chosen near J^; the 

trace-type form is consistent with our gauge in which the cross sections of £ are round spheres relative 

to the unphysical metric. 



D — 2 



Rab — 



(D-1MD-2) 



R~g, 



ah 



D-2 



/i crf (trace) (95) 
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where "trace" denotes the contraction with g cd of the preceding terms. As shown in [22], 
both N a b and 8N a b = ^rN a b(Q) are smooth at J? in our gauge. In fact, eq. (56) from [22], 
which uses the linearized Einstein equations, gives that, at J^ + , 

D — 2 

8N ab = -^nlab + D VL' l h^ b)c n c — ft~ n c n c j ab . (96) 

From this, eq. (189"]) . and using also n c n c = 0(Q 2 ) since u is an affine parameter on scri, 
we obtain 

WV 12 ( 7 ,^f t7 ) = -L f (t a V a u) SN^SN^ + C^-Cih) , (97) 
167r Js„ 



where 

'C(t) 



327T 



and where the natural integration element on C(t) induced by jl a b is understood. 
Next, we evaluate Wjf 12 (7, -^7) . The Raychaudhuri equation on + yields 



-r-tf(A) = -— ^(A) 2 - MAK'(A) - i? ab (A)nV , (99) 
du D — 2 



where, as before, 

= ]^i ah ^ a b (100) 

denotes the expansion of the generators of e^ + , whereas a a b denotes their shear 

^^nfiab = °ab + Jj ~ Hab^ • (101) 

Here, the quantities \i a b and n a = (d/du) a refer to the horizon metric; see ([5]). We take 
a A-derivative of this equation and evaluate at A = 0. Then, since $(0) = = a a b(0) 
for the background metric <? a b(0), and since SR a b = dR a b/d\(0) = by the linearized 
Einstein equations, we have (15$/ du = 0. Since <5-#(0) = on B by our choice of gauge 
(see subsection 2.1), it follows that 5d = on the entire horizon, as we have previously 
claimed. This implies that J*? n (g ab 1ab) = 0. Using this fact, one can show that 

W M2 ( 7 , = ^ l (^Kl Cd ^nlcd - l^n^Klcd) , (102) 

where we have used the axisymmetry (IDT]) of -y a b to replace Jtftlcd by ^fxlcd, where K a is 
the Killing field (fT2"|) normal to the horizon. Using K a = Kun a , we obtain 

W> 12 ( 7 ,^7) = / (^ a V aM ) <Wa cd + £(t 2 ) -B(ti) (103) 



where 



B(f) = -i- / (K a V a u) 7 cd (^n7cd) , (104) 

a-* 71 " J Bit) 
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and where the natural integration element on B(t) induced by fi a b is understood. The 
above calculations motivate the definition of a modified canonical energy £(7, S^(t)) given 
by 

£( 1 ,Y(t)):=£( 1 ,y(t)) (105) 



327T Jq^ 327T j B ( t ) 

Thus, the modified canonical energy £(7, S^(t)) differs from £(7, S^(t)) only by the above 
boundary terms C(t) and B(t). Note that C(t) vanishes when the perturbed Bondi news, 
SN a i,, vanishes. This is seen using eq. ( 196|) . the fact [22] that 7 Q &n a = 0(0,), and that 
n a n a = 0(Q 2 ), which in turn follows since u is by construction an affine parameter in 
Also, B(t) vanishes when the perturbed shear, 5a a b, vanishes, since 5o~ a b = \^ n lah 
on J^ 7 " 1 ", by (jlOip . 5$ = = 1? on Jtf ?+ , and by expression eq. (I149p for 7 ab . In addition, B 
vanishes at B since K°-\b = 0. Since the perturbed Bondi news vanishes as one approaches 
spatial infinit}H it follows that £ (7, J^(t)) -> ^(7, S) as ^(t) -> E. 
The above results establish the following theorem: 

Theorem 1. Let 7 a b be an axisymmetric linearized perturbation that, in addition to the 
properties assumed in the previous section, is asymptotically flat at future null infinity. 
Then the modified canonical energy (11051) has the property that, for t\ < ti 

£( 7 ,^(t 2 ))-£(7,^(ti)) (106) 
= -^T f (t a V a u) 6N cd 6N cd - -L / (K a V a u) 6a cd 6a cd < . 



We remark that Habisohn [20] has considered, in D = 4 dimensions, the second order 
Einstein tensor, and has derived a balance law for the flux through an infinitely extended 
timelike tube approaching null infinity that is similar to that given above in f)106p . al- 
though his balance law does not refer to the flux between fixed times t\ and t 2 as above, 
but to the total flux for all times. 



4 Stability and instability 

As seen is section 2, the canonical energy, S , can be viewed as a quadratic form defined 
on the vector space of smooth, linearized solutions that satisfy our asymptotic conditions 
at spatial infinity and our gauge conditions near the horizon (see subsection 2.1). If £ 
were positive definite on this vector space, it would provide a conserved norm that could 
be used to argue for stability. On the other hand, if £ < for some perturbation, then 
the flux results of the previous section show that the modified canonical energy £ on the 
slice y(t) can only become more negative with time, which suggests instability. However, 
there are a number of obvious difficulties with making such arguments. In particular: 

• As already noted [see (|56|) ]. £ is negative for the "change of mass" perturbation of 
the Schwarzschild black hole, yet Schwarzschild is known to be linearly stable. Such 
"trivial" perturbations must be eliminated from the analysis. 

25 This is manifestly true because we consider perturbations having compact support on E; see foot- 
note mi 
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• £ will, in general, be degenerate on some perturbations, so, even if positive, it cannot 
be expected to be positive definite. Thus, the degeneracies of £ must be carefully 
analyzed. 

In order to properly analyze the degeneracies of £ , it will be very useful to view £ 
as a densely defined quadratic form on a Hilbert space. In the next subsection, we will 
introduce such a Hilbert space V — whose construction and properties may be of some 
interest in their own right — and analyze the degeneracies of £. It should be noted that 
the perturbations in V will have vanishing linearized ADM mass, linear momentum, and 
angular momentum with respect to rotational Killing fields, SM = SPi = SJa = 0, so 
they will automatically eliminate the "trivial" perturbations corresponding to variations 
of the mass or angular momentum in a family of stationary black holes. We will then 
make our stability/instability arguments in subsection 4.2. 

4.1 The Hilbert space V and the degeneracies of £ 

We st arp with the real Hilbert space 

K, = L 2 (E, (T*E) V2 ; hh ) © L 2 (E, (TE) V2 ® h~h ) 

of all square integrable linearized initial data, not necessarily satisfying the constraints. 
Here A 5 is the line bundle over E of densities of weight one half. Thus, the first summand 
denotes the Hilbert space of square integrable symmetric tensors on E, whereas the second 
denotes the Hilbert space of square integrable symmetric tensor densities of weight \ over 
E. In other words, elements of K, consist of pairs (5h ao , 5p ab ) with inner product 



Sih 
Sip 



(tp) ) := L hl 6iKb62hab + hri5 ^ b5 ^° b ■ ( io? ) 



Note that for D = 4,5, the requirement of square integrability will impose faster fall-off 
conditions at infinity on 5h ao (but not 5p ab ) than assumed in section 2 and will exclude 
the possibility of having a nonvanishing 5M, but for D > 6, /C will include all smooth 
initial data satisfying the asymptotic conditions at spatial infinity stated at the beginning 
of subsection 2.1. We may view the symplectic form (j28l as a bounded linear map 
S:K^K: defined h¥\ 



16tt W< s (5 1 g,8 2 c 

where 



Sih 
Sip 



K. 



-h 2 h ar k 



S = \ i c . (109) 



26 Although we will restrict consideration to axisymmetric perturbations when we make our stability 
arguments, the constructions of this subsection do not require axisymmetry, and we will not impose this 
restriction until the end of this subsection. 

27 Note that, in subsec. 12.41 the angles ( | ) as e.g. in (|T2"j) denote the dual pairing between tensors 
and densities, whereas ( | )jc is an inner product between objects of the same density. As a consequence, 
the definition of C* , defined relative to ( | )jc, used from now on, differs from the earlier one in ([T4"f by a 
factor of diag(h~2 ,hz). This also accounts for the difference between S and a in ([77]). 
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Note that S* = —S and S 2 = —J, so, in particular, S is an orthogonal map, S*S = I. 

We would now like to pass to the subspace of /C that satisfies the constraints and 
satisfies our gauge conditions. At first sight, it might appear that this would be a difficult 
task, given that generic elements of /C are not even differentiable. However, the desired 
subspace V can be defined straightforwardly as follows: Let 

W := {^fx9ab £ | X a G C°°,near infinity X a coincides with a rotational Killing 
field plus an asymptotic translation, X a |# tangent to generators of ^ + ](110) 

Here, by Jz?x9ab we mean the corresponding initial data ((J?xg)h, (^x9) P )- The space of 
interest for us is the space V of initial data symplectically orthogonal to W: 



V = W ±s = I {Sh, 5p) G K 



Sh 
5p 



S 



Sh' 
Sp' 



for all (Sh', 5p') G W > . (Ill) 



K 



Thus, V is a closed subspace of /C and, thus, is itself a Hilbert space. Note that since 
S* = —S, we have V ±s = W, where W denotes the closure of W. 

In the following, we will need to work with tensor fields that are in the weighted 
Sobolev spaces = Wp(T). These spaces denote the weakly differentiable tensor fields 
u on E such that p n D( ai . . . D an) u G L 2 for all n < k, where p > is a smooth function 
on S which interpolates between 1 in a neighborhood of B and (x 2 + • • • + x 2 J _ 1 )2 in a 
neighborhood of infinity. The weight factor p n for the n-th derivative is inserted so as to 
force weak derivatives to fall off faster by an appropriate power. The weighted Sobolev 
norm for u G is defined by 



\u\ 



w k 
p 



Y,j^ hl p 2n ( D ^ ■ ■ ■ D ^) u ) Diai ■ ■ ■ Dan)u 



(112) 



We will use the notation C^°(S) to denote tensor fields that are smooth and of compact 
support on E \ B. We will use the notation C°°(S) to denote smooth tensor fields on E 
that can be smoothly extended across B. 
We denote by U the space 



U = P| (V p fc (E) © H/ p fc (S)) 



(113) 



Clearly, U is a subspace of /C, and it follows from the Sobolev embedding theorems that 
U C C°°(^) ©C^E). 

The proof of property (3) of proposition 5 below will require the following lemma: 

Lemma 3. Let Y = (N , Nfi) G © C^°. Then there exists a solution X in the space 
n k (p 2 W^ +1 (i:) ©pWj(E)) C C°°(E) © C°°(E) to the following boundary value problem: 



and X = (N, N a ) satisfies 

N a 



CC*(X) = Y in E, 



iV?7 a , 5$ = 5e = on B = <9E. 



(114) 



(115) 
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Here C and jC* are given by flTOj) and f !74|) and 5$, <5e are the perturbed expansion/volume 
element on B associated with the perturbation (5h,5p) = C*(X), given by 



5ti\ B = h-^(h ab - r) a r] b ) (s Pab - l -5h c c Pah + 2 Pc{a 6h b) c + l -hS J? v 6h ab 

1 



D _ . ^pKbSh c c + p cd 5h cd h ab + pdh ab + 5p c c h ab/ 



B 



(116) 



Se\ B = hh ab -r) a V b )6h ab , (117) 
2 B 

where if is the unit inward normal to B within S. Furthermore, if D > 5, the solution 
is unique, whereas for D — 4, the solution is also unique unless t a is tangent to the 
generators of the horizon (i.e., the black hole is nonrotating), in which case the solution 
is unique up to X a — > X a + ct a . 

A proof of this lemma is given in appendix B. 

The key properties of V are now summarized in the following proposition: 

Proposition 5. 1. Let (Sh,5p) G V. Then (Sh,Sp) is a distributional solution to the 
linearized constraints in the sense that 



C*{X) 



8h 
5p 



= 0, (118) 

K 



for all X a G C °°. 



2. If (5h, 5p) is a smooth element of V, then 5M = 5Ja = SPi = and 8e\s = 5$\b = 
0. Conversely, if (5h, 5p) G K, is a smooth solution of the linearized constraints 
satisfying 5M = 5J A = SPi = and 5e\ B = 5$\ B = 0, then (5h, dp) G V. 

3. The space U fl V is dense in V, where U was defined by (11131) . 

Proof: Property (1) is an immediate consequence of (T75|) together with the fact that for 
any X a G C^°, we clearly have J^xQab G W. Property (2) is an immediate consequence^! 
of lemma [2J 

Our strategy for proving property (3) is to define a bounded projection operator II : 
/C V which takes C^©C^° C L 2 ©L 2 = /C to elements in UnV. Since C£°©C£° is dense 
in /C, it will then follow that U fl V is dense in V. We will produce the desired operator IT 
in two steps. First we define an orthogonal projector II : /C — > Vo onto a closed subspace 
Vo C /C = L 2 © L 2 containing V. The space Vo "imposes the constraints" and the "correct 
boundary conditions" 5e\s = 5$|b = at B, but does not impose 5M = 5. J a = S Pi = 0. 
We will then compose Ilo with a finite co-rank projection operator IIi that commutes with 
n to obtain the desired projection operator EL Since orthogonal projectors and finite 
co-rank projectors are bounded, II is bounded. 

Our prescription for the operator Il is as follows. Let (5h ,Sp ) G © C£°. Then, 
obviously, Y := C(5h ,5p ) is smooth and of compact support. Let X be the solution to 

CC*(X) = Y (119) 



28 0ur asymptotic conditions at infinity are now that (Sh, 5p) £ K. rather than the conditions that 
follow from those stated in the second paragraph of subsection 2.1. It is easily verified that the proof of 
lemma [5] continues to hold with our present asymptotic conditions. 
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given by lemma [3] above. We defined 

TT ( 5k ° 
110 Uo 



5 h 



- C*(X) . 



(120) 



It follows immediately that Uo(5h , $Po) £ that it satisfies the linearized constraints, 
and that it satisfies 5$\b = = 5e\s [cf. eqs. (I116p and fl!17p ] . Since we have 



\C*(X)\\% = (X\Y) 



X 



c 



5 h 
$Po 



C*{X) 



5 h 
Sp 



< \\C*(X)\\ K 



5 h 
$Po 



;m) 



(where eq. fll64p of Appendix B was used in the first equality) it follows that n is bounded 
and, hence, its action can be extended to /C. We now show that Il = ITq and ITq = n o , 



Il , we let #o 



so Ilo is an orthogonal projection. To show that IIq 

:= (5h' ,5p'o) be smooth and of compact support. Then 

<ft J %\%) K = (* W K -(C*{X)\* ) K 
= <¥ l*o>,c - ( x \ Y ')>c = (%\%) K ~ (C*(X)\C*(X')) K 



(5 ho, 5po) and 



;i22) 



where X, X' ,Y,Y' are defined as above and, in the last step, we again used (11641) . The 
expression on the right is manifestly symmetric, thus proving the claim. Similarly, we 
have 



[U %\U %) K 



(123) 



thus showing that IT5 = n o . 

Let Vo denote the closure in K of the image of ©C^° under n . We now show that 
V C V . Suppose that (5h, Sp) G V is orthogonal to U (C^° © C£°). Then inserting (11201 . 
we see that 







5h 
Sp 



8 h 
$Po 



K 




C*(X) 



(124) 



K 



The second term vanishes, because X G p 2 W 2 © pW^, in which C^°(S) © C^°(S) are 
dense. Then, from the boundary conditions of X (see lemma [3]), X can be approximated 
by X n G W such that C*X n C*X in L 2 © L 2 . Hence, since (6h,Sp) in W Ls , it 
follows that ((6h,Sp)\ C*{X))x, = 0. Thus, we see that = ((5h, 6p)\(6ho, 5po)}ic for any 
compactly supported (Sh ,Spo), which clearly means that (Sh,Sp) = 0. Thus, the closure 
in K of n (C °° © C °°) contains V. 

On the other hand, since all elements (Sh, dp) G n (C^° © C^°) are in U (and, hence, 
are smooth), satisfy the linearized constraints, and satisfy 5$\b = = 5e\s, it follows from 
property (2) of this proposition that (8h, 8p) G V if and only if 5M = 5Ja = 8Pi = 0. 
Therefore, to obtain the desired projection map we need only compose Ilo with a finite 
co-rank projector Hi defined as follows: Let / = 1, . . . , k denote a finite collection of 



29 See footnote [STJ 
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smooth vector fields that vanish in a neighborhood of B and coincide, respectively, with t a , 
ipA a , and (d/dx l ) a in a neighborhood of infinity. Let := (Sih, dip), i — 1, . . . , k be an 
(arbitrarily chosen) collection of perturbations in no(C^ ©C^°) for which SH^fifj) = dij. 
The desired projector ITx is defined by 



7=1 





(125) 



and the desired projector II onto V is then defined by II = Ilillo. By our arguments above, 
the image of U.(Cq° © C^°) under II is contained in U and is dense in V thus proving that 
U fl V is dense in V. □ 

Remark 1: It is possible to strengthen property (3) using gluing techniques to show 
that even (C °°(S) © C£°(£)) n V is dense in V, where by definition Cg°(£) © C °°(E) C U 
consists of all smooth initial data that vanishes in a neighborhood of spatial infinity. A 
statement and proof of this strengthened result is given in Appendix C. 

Remark 2: All elements of W are smooth solutions to the linearized constraints and 
satisfy 5M = 5Ja = SPi = and 8e\s = 8$\b = 0. Thus, by property (2) of the above 
Proposition, we have W C V and, hence, W C V. On the other hand, since W and V 
are symplectic complements of each other in /C, it follows immediately that an element 
7 = (Sh, dp) G V is such that Ws(g; 7', 7) = for all 7' G V if and only if 7 G W. Using 
the inequality f)163p . it can be shown that smooth elements of W must lie in W. Thus, 
for any 7 G W fl V, we have W-s(g; 7', 7) = for all 7' G V if and only if 7 G W. 

We now restrict consideration to axisymmetric perturbations, =^ A 7 = (see flo^]) ). We 
can redefine the spaces /C, W, and V with the word "axisymmetric" suitably inserted, and 
all of the results of this subsection and their proofs continue to hold without modification. 
In order not to make our notation more cumbersome than necessary, we shall continue to 
use /C, W, and V to denote the axisymmetric versions of these spaces. Thus, even when 
not stated explicitly, axisymmetry should be understood in all statements below. 

We now view the canonical energy S, defined in section 2, as a quadratic form on the 
(axisymmetric) Hilbert space V, with dense domain T = U H V. We are interested in 
finding all of the elements 7 G T on which £ is degenerate, i.e., for which 

£( 7 ' )T ) = for all 7' GT (126) 

Since £(7', 7) = W^{g\ 7', ^7), it follows immediately from Remark 2 above that £ is 
degenerate on 7 G T if and only if ££ t "j G W. But, in the presence of axisymmetry, 
this is precisely the condition that 7 is a perturbation towards a stationary black hole 
(see Definition 2.1 at the end of subsection 2.3). Let T' denote the space of equivalence 
classes of elements of T, where two elements are equivalent if and only if they differ 
by a perturbation towards a stationary black hole. We have proven the following: The 
canonical energy £ is well defined, non- degenerate quadratic form on T' X T' . 



4.2 Stability and instability arguments 

For axisymmetric perturbations of a stationary black hole, we have just shown that £ is 
well defined as a quadratic form on the space T' of smooth solutions satisfying our gauge 
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conditions with 5M = 5Ja = 5Pi = modulo perturbations towards a stationary black 
hole. Since £ is non-degenerate on T', only the following two cases can occur: 

• Case (a): £ is positive semi-definite on T and hence is positive definite on T' . 

• Case (b): There exists 7 G T such that £(7,7) < 0. 

We now analyze the behavior of perturbations in these two cases. 

Case (a): £ provides a positive definite conserved norm on T' . This precludes existence 
of "growing modes" in T' . But elements of T can be expressed as the sum of a represen- 
tative element of T' and a perturbation towards a stationary black hole. Since the latter 
perturbations are manifestly stable, this shows stability for perturbations in T. However, 
by Proposition 5, T is dense in V. Now, a general axisymmetric solution to the con- 
straints in K, can be written as a sum of an element of V and a finite linear combination of 
representative "change of mass," "change of angular momentum," and "change of linear 
momentum" perturbations. The desired "change of linear momentum" perturbations can 
be chosen to be ££y t g where Yi is an asymptotic boost. These are manifestly stable. If the 
black hole is part of a family parametrized^! by (M, J a), then the "change of mass" and 
"change of angular momentum" perturbations can be chosen to be the perturbations to- 
wards other members of this family. These are also manifestly stable. Thus, we conclude 
that the black hole is stable to a dense set of axisymmetric perturbations in the subspace 
of /C comprised by solutions to the constraints. 

Furthermore, the results of section 3 strongly suggest that for perturbations 7 G T, £ 
and £ should decay to zero on the slices y{t) as t — > 00, and thus, at late times, 7 should 
approach a perturbation towards a stationary black hole. 

Case (b): Let 7 G T be such that £(7,7) < 0. We obtain a contradiction with the 
possibility that 7 approaches a perturbation, 70, towards a stationary black hole on <9*{t) 
as t 00 as follows. If 7 approached a stationary solution, then the Bondi news and 
the shear of the horizon should approach zero at asymptotically late times, in which case 
as t 00 we should have £(7, S?(t)) —> £(7, J^(t)) —> £(70) = 0. However, this is a 
contradiction because for all t we have £(7, y{t)) < £(7, S) = £(7, S) < 0. 

Although this contradiction does not prove instability, the fact that \£('~f, y(t))\ can 
only increase with time suggests that the amplitude of the perturbation does not decrease 
with time. But then the fluxes through the horizon and null- infinity should 

also not decrease, causing (£(7, y{t))\ to increase further, etc. It therefore seems highly 
plausible that £(7, y(t)) — > —00 as t — > 00 and that the amplitude of the perturbation 
grows without bound, as has been previously argued in a different context in [HI [T2J dSl 

nn. 

Thus, we conclude that in Case (a) the black hole is stable, whereas in Case (b), the 
black hole is unstable. In other words, we have argued that the necessary and sufficient 
condition for black hole stability is that the canonical energy £ be positive semi-definite 
on the space T = U fl V of smooth solutions with 5M = 5Ja = SPi = that satisfy our 
horizon boundary conditions 5$\b = o~c\b = 0. 

30 Even if the black hole is not part of a family parameterized by (M, J4), it should be possible to show 
existence of a basis of stable "change of mass" and "change of angular momentum" perturbations. 
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We now briefly indicate some possible ways of improving the above arguments. To 
prove stability in Case (a) , we would like to establish boundedness of 7 itself rather than 
have boundedness/conservation of a norm on 7. However, this is clearly not possible 
until we have imposed suitable gauge conditions on 7 that uniquely fix its evolution. 
Indeed, our present gauge conditions restrict 7 only near infinity and near M )+ ^ so they 
allow pure gauge perturbations for which 7 becomes arbitrarily large. If suitable gauge 
conditions are imposed on 7, then it is not implausible that £ would provide a conserved 
norm on 7 G T that would be equivalent to a Sobolev norm. The canonical energy of 
the perturbation (J£ t ) kr y would similarly provide a family of conserved norms on 7 G T 
that plausibly would be equivalent to higher Sobolev norms. If so, boundedness of 7 itself 
would be proven. 

It does not seem feasible to try to directly convert the above argument for instability 
in Case (b) into a mathematically rigorous proof. However, it is possible that a proof 
of instability in Case (b) [as well as stability in Case (a)] could be obtained along the 
following lines. Again, to get started, it should be necessary to suitably fix a gauge so 
as to obtain deterministic dynamics. Consider, first, the case of perturbations off of a 
static black hole. On account of the time reflection symmetry of the background, £ will 
be invariant under time reflections, so the integral expression for £ in terms of initial data 
will not contain any "cross-terms" between Sp and 8h. Thus, £ can be written as a sum of 
a "kinetic energy" (quadratic in 8p) and a "potential energy" (quadratic in 8h). It is likely 
that the kinetic energy is always positive-definite. For the case of spherically symmetric 
perturbations (with matter fields) off of a static, spherically symmetric background, it was 
shown in [35] that the kinetic energy could be used to define an inner product such that 
the potential energy could be expressed as the expectation value of a self-adjoint operator 
A that appears in the dynamical evolution equations. In this case, it can be seen directly 
from the dynamical evolution equations that positivity of A is equivalent to stability. 
However, since the kinetic energy is positive, positivity of A is equivalent to positivity 
of £, thus rigorously establishing the equivalence of dynamical stability and positivity of 
£. It is not inconceivable that similar results could be proven for perturbations of an 
arbitrary static black hole, although, obviously, a number of technical details would have 
to be sorted out. A similar strategy could also be applied to the stationary- axisymmetric 
case — making use of the t — ipA. reflection isometry to decompose initial data into its "time 
symmetric" and "time antisymmetric" parts — but it would be less obvious in this case 
that the "kinetic energy" would have to be positive. 

5 Proof of the Gubser-Mitra conjecture 

Up to this point, we have restricted our considerations to the stability of black holes. In 
this section, we extend our considerations to black branes of the form ([1]). By making 
the same type of dynamical instability argument as given in subsection 4.2 above, we will 
show that if a family of black holes is thermodynamically unstable, then the corresponding 
family of black branes must be dynamically unstable. 

To begin, suppose we have a family of stationary, axisymmetric black holes, g a b(M, J a), 
in D > 4 dimensions, that satisfy the asymptotic and gauge conditions of subsection 2.1. 
Let £ = (Mq, Jqa) denote the parameter values of a particular black hole in this family, 
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and consider the one-parameter subfamily 



9ab(X) :=g ab ^ + \v) , (127) 

where v is an arbitrary vector in the parameter space (M, J a)- It is obvious that for this 
family, we have 

d 2 M . , d 2 J A , . 

^(0) = ^(0)=0. (128) 
Thus, by fl53|) . for the perturbation determined by v we have 



On the other hand, we have 



d 2 A 



(0) = H.ess A \dv,v) , (130) 



dX2 

where Hess^i* denotes the Hessian of A [see eq.([2D] at parameter value £. Thus, we see 
that one can find a perturbation to a black hole in the family that makes £ negative if 
and only if the Hessian Hess a has a positive eigenvalue. 

However, a nontrivial perturbation to a black hole in the family clearly does not have 
5M = 5Ja = 0, so a negative value of £ — or, equivalently, a positive eigenvalue of Hess a— 
does not provide any information about stability. However, consider now the (D + p)- 
dimensional spacetime ^# = ^ x T p with metric g a b defined as above in eq. (JTJ) with each 
Zi a 27r /-periodic coordinate parameterizing the corresponding "extra-dimension". (Here 
and in the following, a tilda dentotes a quantity associated with the (D + p)-dimensional 
spacetime This metric represents a "uniform black brane", where uniform refers to 
the fact that each (d/dzi) a is a Killing field of this spacetime, and where brane refers 
to the fact that the horizon cross section is now B = B x TP. A black brane spacetime 
is "asymptotically Kaluza-Klein (KK)" rather than asymptotically flat, but the ADM 
conserved quantities (M, P,J,C) can be defined as in the asymptotically flat case (see 
eq. (f3"7]) ). In addition, we have conserved quantities ("KK charges"), T, associated with 
the asymptotic symmetries (d/dzi) a , given by 

5fi-.= [5Q 9/dz i(g) - i 9 /dz^(g,Sg)] . (131) 

J oo 

We now prove the following proposition: 

Proposition 6. Let g a b(M, J a) be a family of black holes that is thermodynamically 
unstable at (M , Joa), i-e., there exists a perturbation v within the black hole family for 
which £ < 0, which will be the case if Hess^ has a positive eigenvalue at (M , Joa)- Then, 
for any black brane corresponding to g a b(M , J a) via eq.flTJ with sufficiently large I, one 
can find a perturbation for which £ < and 5M = 5Ja = 6 Pi = Ti = 0, = 5e\^ = 0. 



Proof: For notational simplicity, we prove the proposition for p = 1 extra dimension, 
the general case is completely analogous. Let (5p ab , 5h a b) denote the initial data of a 
perturbation v in the black hole family for which £ < 0. Without loss of generality, 
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we may assume that our Cauchy surface E is chosen to be maximal in the background 
and perturbed spacetime^]], so that p ab h ab = 0, and 5p ab h a b + p ab 5h a b = 0. We wish to 
construct initial data (5p ab , 5h a b) in the (D + l)-dimensional black brane spacetime with 
£ < for which 5M = 5J A = 5P i = 5T = 5A = and for which 5$ = at B. 
Our ansatz for (5p ab , Sh ab ) is 

Sp ab =[5p ab + ^yP a V + h'HD^XV + £>(°C(£) 6) - ^(h ab + (jLy(jL) b )D c X c 

+ ik/l (£ + ^ iA;/Z C (£) (a (£) b) - ^ ik/l Ch ab }) e^ 1 , (132) 



where k G Z and 0, £, X a are tensor fields on E (so they are independent of z, and 
tangent to E). The terms in 8p ab involving £, X a may be written alternatively as D^ a X b ^ — 



±h ab D c X c , where 



X a = (C (|-) a + X a ) e* 2/i . (133) 



If G Z is not zero, then it is easily seen that the surface integrals defining SM, 
5 J a, 5Pi, 5T and 5 A vanish because the z- dependence is e lz ' k ^ 1 . 

The tensor fields <ft, (,X a in our ansatz f l X 3 2 j) are now chosen so that the linearized 
momentum and Hamiltonian constraints are satisfied. It can be seen that the satisfaction 
of the linearized momentum constraint is equivalent to 

D^D^XV - ^h ab D c X c ) = Uk/l h-hh cd p cd e lk z / l (£) b . (134) 

This is an equation on E, but because all tensor fields have the same dependence e tk ' z ^ L 
on we can rewrite it as the following system of equations on E: 



D a (p( a XV - ^h ab D c X c ^j - i k 2 /l 2 X b + ik/l D b ( = 0, 

D a D a ( - 2^1 e/f C + £j_l ik /i D - Xa = ik/l h-hh cd p cd 



(135) 



The satisfaction of the Hamiltonian constraint is equivalent to the linearized Lichnerowicz 
equation for <fi, 



-D a D a + h- L p abP ab + k l /V 



(D-2) 



(136) 



Together, the eqs. (I136p . fll35p form a system of linear, inhomogeneous, elliptico PDE's 
for the unknown tensor fields <f), (,X a on E. We need to impose boundary conditions 



31 Existence of a maximal slice in the background spacetime is shown in (3J. For the perturbed spacetime, 
existence of a slice that preserves maximality follows from the fact that the perturbed lapse required to 
achieve this condition satisfies a suitable elliptic equation. 

32 Indeed, the principal symbol of the equation for X a in (JT35)) is a(£) a b = h cd ^ c S a b - (1/D) h ac ^ b , 
which is an invcrtible linear map in each tangent space TE for any T*£ 3 £ a =/= 0. 
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at B. To find these, we consider the linearized expansion 5d on B in the black brane 
spacetime. We would like this to vanish. It is calculated using eq. ( I116P as 

<W =( - ^ h-hab(h ab - r/V) - rfD a <l> + ~ cf>JZ> v h ab (h ab - V a V b ) 

- rfrfD a X h + 1 D C X C + 1 ifc/Z c) e*^' 

- TfA^ + 2 #0 - 77 a 77 6 D a X 6 + 1 D C X C + 1 ifc/Z () e**/ 1 (137) 

- ^,0 + 2 #0 - ^°X ) + 
■i X a ) + 1 X? a [(/i" 6 - tjV)^] + 1 ifc/i C) e**-*/' at B, 



where ?? a is the normal to B within S, D a is the derivative operator intrinsic to B, and if 
is the trace of the extrinsic curvature of B within S. In the first line, use has been made 
of the fact that p a a = and 5$ = on B in the original spacetime. In the second line, 
we used that = i? = {h ab — T] a r] b )(h~2p ab + ^<5f v h ab ) on B in the original spacetime. In 
the last step we have, without loss of generality, extended rf of B such that it is geodesic. 
Because B is the bifurcation surface of a Killing horizon, H = on B. In view of this, 
5$ = on B will follow if we impose on 0, (, X a the boundary conditions 

^0|b = , 

(h ab - 7] a 7] b )X b \ B = , , 
v / / 138 

^(X%)| B = 0, 

C|b = o. 

In other words, we impose Neumann conditions on 0, Dirichlet conditions on (, Neumann 
conditions on the orthogonal components r] c X c of X a , and Dirichlet conditions on the 
tangent components (h ab — r] a r] b )Xb of X a at B. We need to analyze the existence and 
properties of the solutions to (I136p . (11351) . The existence of a solution ((,X a ) £ pW^(S) 
(and in fact £ /?W* for any k) follows from lemma H] in appendix B, because (11351) is 
equivalent to (11341) . and the lemma gives a smooth solution X a £ pW^(E) to the latter 
equation. Next, we analyze the existence and properties of the solutions to (I136p . We 
need to distinguish the cases D = 4, 5 and D > 6. 

D > 6: In this case, the right side of H136[) is in L 2 . Existence of a smooth solution 
in the Sobolev space pW^ (and in fact £ pWff any fc) satisfying Neumann boundary 
conditions can be proven as usual by considering the weak formulation of the boundary 
value problem 



(D-2) 



D a <f)D a ^ + h- x V ab Vah ^ + k2 / [2 ^) hh 

J (k 2 /l 2 hhh a a ^-2p ab D a X b ^j 



(139) 



for all smooth %p with compact support in E, and making use of the Poincare-type in- 
equality 

c / p- 2 ip 2 hh < / D a ^D a ip hh , (140) 
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for some c > 0. Since these arguments are standard and very similar to those in the proof 
of lemma [3] in appendix B, we do not repeat them here. 

It remains to be shown that the normalized canonical energy £ = 1/(2ttI) W^(5g, J£t&g) 
on the (D + l)-dimensional black brane spacetime approaches the canonical energy £ < 
of the original black hole perturbation for / — > oo. For this, it is important to understand 
how the norms of (pi, Q, Xf 1 behave for large / (here we put an "/" on the tensor fields to 
indicate their dependence on /). We use that the source term in (I134p is of order 0(p D ^ 1 ), 
so p times the source is in WpiT,) in D > 6 and tends to zero in that normal From 
lemma H] of appendix B we infer that the L 2 (S)-norms of p~ x Xi a , p~ 1 DbXi a , pD a Df,Xi c 
tend to zero of order 0(1//) as / — > oo. This translates into the statement that the 
L 2 (S) norms of p~\ h p^Xf, D b Xf, D a Q h p/l D a £ h p/l D a X t b , l/l 0,1// as well as 
pD a D b (i, pD a D h Xf tend to zero of order 0(1/1) as / — > oo. It follows immediately 
from f ll36p that the L 2 -norm of l/l (pi remains bounded as / — > oo. Then, combining f)139p 
with (I140p for = ijj — (pi, we find that p~ x (pi, D a (pi — >■ in the sense of L 2 as / — > oo. 
By standard bootstrap and scaling arguments for elliptic PDE's of the nature under 
consideration (see e.g. [6] Appendix A for a discussion and references), it follows similarly 
P D a D b (Pi -> in L 2 . 

Using now the conditions dp^hab+p^Shab = 0, p c c = (from which follow the same for 
the corresponding tilde tensor fields on S), the boundary conditions at B, and eliminating 

33 In D = 4, 5 this is still true if we impose a gauge, as below, that 5h a b — jj^ihabSh is of order 
0{p- {D - 2) ), because p a a = 0. 
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terms that are total derivatives with respect to z, we obtairo from (l86j) 
16/T £ — 167r £ = 

Re / N^R ab (h)6 1 h c c *S 2 h ab + ^R ab (h)6 1 h ab *6 2 h c c - 2R ac {h)8 l h ab * 5 2 h c b - 

\b x h ac *D a DJ> 2 h b b - -D a D c b x h b b *b 2 h ac - U x h ac *D b D b b 2 h ac - ^5 2 h ac *D h D b 5 x h^+ 
X -hh ac *D h D a b 2 h c h + X -D b D a 8 x h b *b 2 h ac - ^DaD^h^Sihc)- 
^D a {S l h' M D b S2hc c + 5 2 h a *D b 8 l h c c *) + ^D a D a (8 1 h b b *5 2 h c c )+ 
D a (8ih c a *D b 5 2 h cb + S 2 h a c D h S x h A *) + l -D a {8 x h b *D h 5 2 h ac + SJifD^h"*)- 
^D a (5 1 h/*D b 5 2 h ab + 5 2 h d d D b 5 1 h ab *)} h?+ 

Re J N{±p ab Pab 6 1 h a a *6 2 h b b - ^p c bP a \5 1 h d d *5 2 h ac + S 2 h d %h ac *)+ 
26 lP ab *6 2Pab + 2 P ac p b %h ab *6 2 h cd - ^5thi*5 2 p ab p ab - 
p 2 h d %p ab * Pab + Ap c b {5 l h a *8 2 p oh + 5 2 h ac 5 lP ab *)} hr^- 

Re J iV a | - 6 lP bc *D a 5 2 h bc - 5 2 p bc D a 5 1 h bc * + 25 lP cb * D b 5 2 h ac + 
25 2 p cb D b 5 x h a * + 8 1 h ac *D b 5 2 p bc + 6 2 h ac D b 5 lP cb *- 

V cb {b x h ac *D b b 2 h d + 6 2 h ad D b 5 1 h d *) + ^p cb (6 1 h ad *D%h cb + feWA*)}- 

J B 

Re / iv{(D C)*D a C + -^(D a X a YD b X b + 2{D ~ 2 1} ' k 2 /l 2 (*( + k 2 /l 2 X* a X a + 

X - k 2 /l 2 (§h ab - ^—th^* (5h ab - -^—^<t>h ah )} + 
Re J N a [ik/l {5h ac -^—^(j>h ac ){DX + ik/l X c )*} hh . 

(141) 



34 Note that £ below is extended to complex-valued perturbations in such a way that it is anti-linear 
in the first argument; of course we could equally work with real valued perturbations by taking the real 
and imaginary parts of eqs. (|132p . 
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where, in this equation, * denotes complex conjugation. Here we have defined 



5iK b = - — (ph ab 

6 2 h ab = 25h ab - jy—^4>hab 
1 



(142) 



S lP ab = —±— p ab <f) + h^{D {a X b) -^(ik/lC + D C X C )} 

S 2 p ab = 25p ab + jy—[P ab( P + h^{D {a X b) - —(ik/l ( + D C X C )} . 

From this expression, together with the decay of the norms of the tensor fields Q,(pi, Xf 
which we have described, it follows^] that the terms on the right side of (11411) go to zero 
as I — > oo. Thus, we have £ — > £ < as I — > oo, as we desired to show. 

D = 4,5: In this case, the right side of HI 36[) need not be in L 2 , and the previous argument 
for obtaining a solution <p does not work directly. But in this case, we can apply the 
following modification to the above argument. We now write (p = xbh a a + v, where x is 
equal to 1 near infinity, and equal to in an open neighborhood of B. Imposing 



on <p gives an equation for v of the same type, but with a source in L 2 (and with the 
boundary condition J£ v v — at B). Existence of a smooth solution v G W 1 again follows, 
and therefore existence of a solution <p follows, too. Now let I — > oo and denote the 
corresponding sequence of solutions by vi. From the equation satisfied by vi it is easy to 
see that D a vi is Cauchy in L 2 , and, by (I140p . it follows that p~ x v\ also is Cauchy. Hence 
p~ x vi,D a v\ are convergent sequences in L 2 , with limit p~ x v G W 1 . But it is easy to see 
that v must be equal to —xbh a a since they both satisfy the same boundary value problem, 
so (pi again satisfies p~ l (j>u D a (f>i, pD a D b (f)i — )■ in L 2 . It then follows that all terms (I14ip 
again go to zero. The only terms for which the argument is somewhat different than for 
the case of D > 6 are the last two involving 5h ab — -jjzi§ih ab , because 5h ab is not in L 2 
now. For that term, we write 5h ab — j^zi<t>ih ab as (6h ab — j^Y5h c c h ab ) + j^zi(Sh c c — (p^h^ 
Although 5h ab is not in L 2 in dimensions D = 4, 5, we may assume that we are in a gauge 36 ! 
such that 5h ab — jjzi5h c c h ab falls off faster by one power of p" 1 , implying that it is in 

L 2 . Furthermore, we have (-D a D a + hr x VabV ab )<pi = \jjE§j{k 2 /l 2 ((pi - 5h a a ) - 2p ab D b Xf]. 
Because p~ x (pi, D a (pi, D b Xf — > in L 2 , it then follows (multiply by ((pi — Sh a a )* and 
integrate by parts) that l~ x ((pi — 5h a a ) — > in L 2 . Thus, again £ — > £ < 0. 

Since 5A = 0, we can further change the initial data by a gauge transformation 
[i.e. adding aC*(Y) with a suitable vector field Y a tangent to B, see eq. (177)) ] so that 
5e\g = 0. For example, take Y a = j£-r (ik/l)^ 1 <P(-§^) a e tk ' z ^ 1 on S near B, and vanishing 

in a neighborhood of infinity, and extend it off S such that JzffF a = 0. By proposition 3, 
this will not change £. □ 



35 It is essential here to use the conditions p a f, = 0(p~( D ~ 2 ' 1 ), R a b(h) — 0(p^^ D ^ 1 ' > ) as well as Sh a b = 
0(p~ l > D ~ 3 ^), 5p ab = 0(p~ l > D ~ 2 ^). One can see that (/>;,^,X ; a and their derivatives only enter in the 
combinations such that their norm tends to zero in L 2 , or such that the additional 1/Z-factors give 
convergence to 0. 

36 For example, for perturbations from a Schwarzschild to another Schwarzschild black hole, putting the 
metric perturbation into "isotropic coordinates" near infinity would provide such a gauge, and similarly 
for perturbations within the Kerr family. 
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Essentially all of the analysis of sections 2, 3, and 4 can be applied straightforwardly 
to black branes. One important modification, however, is that in the definition of "per- 
turbations toward stationary black branes," allowance must be made for the possibility 
that the horizon Killing field K a (X) may correspond to a linear combination of the Killing 
fields (d/dz t ) a in addition to t a and Thus, definition 2.1 must be modified by in- 
cluding linear combinations of the Killing fields (d/dz l ) a in eq. (l64"|) . Correspondingly, in 
Proposition 4, the requirement that 71 also satisfy STi = must be added to the hy- 
pothesis. In section 4, the definition of W, eq. flllOp . (and, hence, V) must be similarly 
modified to allow X a to coincide near infinity with a combination of the Killing fields 
(d/dz l ) a as well as rotational Killing fields and spatial translations. Correspondingly, in 
statement (2) of Proposition 5, the condition STi = must be added to the conditions 
5M = 5Ja = 5 Pi = 0. Aside from these changes, the results of the previous sections 
apply. 

We have just shown that for a black hole family such that the Hessian matrix Hess^ 
has a positive eigenvalue, there exists a perturbation of the corresponding black brane 
that satisfies 5M = 5Ja = SPj = 5% = together with 5$\g = 5e\^ = 0, and, for 
sufficiently large I, has £ < 0. Since the flux analysis of section 3 for asymptotically flat 
spacetimes carries over to the asymptotically Kaluza-Klein case, we may make the same 
arguments as given for Case (b) of section 4.2 to conclude that the black brane must 
be unstable. Indeed for D > 6 the perturbation we have constructed lies in the black 
brane analog of V, so we may make precisely the same arguments for instability. However, 
these arguments can also be made for the case D = 4, 5, since the black brane analog of 
Proposition 4 of section 2.3 yields a contradiction with the perturbation approaching a 
perturbation towards a stationary black brane at late times. Thus, we conclude that for 
any family of black holes that are thermodynamically unstable, the corresponding family 
of black branes is dynamically unstable to sufficiently long wavelength perturbations. 

6 Equivalence to the local Penrose inequality 

In this section, we show that for black holes, the satisfaction of the local Penrose inequality 
is equivalent to the positivity of canonical energy for perturbations with SM = 5Ja = 
SPi = and hence, by our previous arguments, is equivalent to dynamical stability. 

Let g a b{M, J a) be a family of stationary, axisymmetric, and asymptotically flat black 
hole metrics on M. Let 5w>(A) be a one-parameter family of axisymmetric metrics such 
that g a b(0) = g a b(M , J a). Without loss of generality, we may assume (by applying 
an asymptotic Lorentz boost, if necessary) that the linear momentum of g a b{X) vanishes, 
•Pi (A) = 0. Let M(A), Ja(A) denote the mass and angular momenta of g a b(X). Let g a b(X) = 
g a b(M(X), Ja{X)) denote the one-parameter family of stationary black holes with the same 
mass and angular momenta as g a b(X). Let -4(A) denote the area of the apparent horizon 
of (?ab(A) on an initial slice and let -4(A) denote the area of the event horizon of g a &(A). 
Then the following proposition holds: 

Proposition 7. There exists a one-parameter family g a b{^) for which 

„4(A) > A{X) (143) 
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to second order in A if and only if there exists a perturbation j' ab of g a b(M , J a) with 
SM = 5J A = 5Pi = such that 8{j') < 0. 



Proof: On account of our gauge condition 5$\b = 0, the apparent horizon of g a bW 
coincides with B to first order in A. To find the apparent horizon of g a &(A) to second order 
in A, we would, in general, need to displace B by a second order gauge transformation 
similar to (jSJ). However, since the extrinsic curvature of B vanishes for the background 
metric <?o&(0), the area will be unchanged to second order under such a displacement. 
Consequently, to second order in A, we have -4(A) = 4(A), where 4(A) denotes the area 
of B in the metric g a b(\)- Similarly, to second order in A, we have .4(A) = 4(A). 

Obviously, we have 4(0) = 4(0), as the two families coincide for A = 0. Likewise, by 
the first law of black hole mechanics, we have 



dA, , 8tt 



dM 
~d\ 



dA, , 

« (0) 



(144) 



Thus, what matters is the second derivative of the areas. Let 7 a & = (dg a b/d\)\\=o and let 
lab = (dg a b/dA)\\=o. Since the families have the same mass and angular momenta, our 
second variation formula f )53|) gives 



K 

8^ 



d 2 A 
~d>? 



(0) 



d 2 A 
~d>? 



(0) 



5(7,7) -5(7,7) 

-5( 7 , ,7') + 25( 7 , ,7) • 



(145) 



where we have defined Y ab = 7afe — lab- However, j' ab is clearly a linearized perturbation 
with 5M = 5Ja = SPi = 0, and 7a& is manifestly a perturbation towards a stationary 
black hole. Therefore, the second term vanishes by Proposition 4 of section 2.3. Thus, we 
have 

d , . d A , . 

-(0) - -ttt(o) 



K 

8tt 



d\ 2 



d\ 2 



-5(7', 7') 



(146) 
□ 

is incompat- 



from which the proposition follows immediately. 

As explained in the introduction, the satisfaction of the inequality 
ible with the stability of the family of black holes g ab (M, J a). It is therefore reassuring 
that we have found that satisfaction of this inequality is equivalent to the condition that 
we have obtained in this paper for dynamical instability. 
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A Proof of Lemma 



2 



Lemma 2. Let 5g ab be a solution to the linearized Einstein equations satisfying our 
asymptotic flatness conditions and our gauge conditions (J7J) and (ITTj) at B. Suppose in 
addition that 5 A = (so that, by (ITT]) , we have 5e|# = 0) and that SH X = for some 
asymptotic symmetry X a . Then Ws(g; 5g,Jf^g) = for all smooth £ a such that (i) £ a |# 
is tangent to the generators of Jrff + and (ii) £ a approaches a multiple of X a as p — > oo. 
Conversely, if 5g ab is smooth and asymptotically flat and if Ws(g; Sg, Jzf^g) = for all 
such £ a , then 5g ab is a solution to the linearized Einstein equation with 8$\b = Se\s = 
at B and with 5Hx = 0. 

Proof: For a solution 6g ab , our fundamental variation formula ( 1551) yields 

Wx{g;6g,JZ>tg) = [ [5Q^g) - i^(g;5g)] - I [5Q^g) - k6(g;6g)} 

Joo J B 

= - [ [SQtbj) - n9(g; 6g)\ + SH^ . (147) 

J B 

Since £ a — > cX a for some constant c, we have 5H^ = c5Hx- 

We now evaluate the boundary term from B. To do so, we write n a = (d/du) a ,l a = 
(d/dr) a , where (u, r) are the Gaussian normal coordinates of ([5]). The condition that C, a \B 
is tangent to the generators of implies that we can decompose it as 

£ a = fn a + uZ a + rY a , (148) 

where Y a ,Z a are smooth but otherwise arbitrary. In Gaussian null coordinates the 
metric perturbation takes the form 

8g ab = -2 V (a u (r 2 5aV b) u + r 5/3 b) ) + 5p ab (149) 

where n a 5j3 a = l a S/3 a = n a 8p ab = l a 5fi a b = 0. 

We now calculate i^9(g; Sg)\B [see eq. (ITT]) ]. For the pull-back to B, we have 

{k9) ai ...a D _ 2 = f £ ai ...a D _ 2 , (150) 

where e ai _ aD 2 is the unperturbed intrinsic volume form on B. From the definition, 
eq. ( Tl8]) . of v a we have on B 



v a n a = g fh n e V f 5g he -g he nfV f 5g he (151) 
= - g f h n e V f (2rV (h u5/3 e) ) + <f VV /(2rV \ h u50 e) ) + /VV/fy a6 - g eh n^ f 5p ab 

The first two terms on the right hand side can give a non-zero result on B (i.e. r = 0) 
only if the derivative V/ acts on r, but in that case the tensor contractions are seen to 
give a zero result. For the last two terms on the right side, we get, again on B, 

v a n a = 2n if l h) n e V f 6p he - 2nH h) n f V f 5p he + /j fh n e V f 5p he - //VV f 5/x he 
= -fi eh n^V f 5p he -fi fh (V f n e )5p he 

= -p ab J?Jii ab + l -5p ab J£ n p ab 

= -p ab ^Jp ab = -26$ . (152) 
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Here, we have used in the first line the expression for g ab in gaussian null coordinates (0), 
whereas to go to the second line we used that at least one vector n a , l a gets contracted 
into 5fi a b when all vectors are pulled through the derivative, and this gives a zero result. 
To go to the third line we used standard expressions for the Lie- derivative of a tensor, 
and that u J a c u J b d V c n d = ^J?? n fi ab on B, which follows from fl5J). To go to the fourth line we 
used that J£ n Hab = 0, as the background is stationary. Thus, we have shown that 



ud = -- 



B 



in 



f8#6 



B 



Next, we calculate 5Q^. When pulled back to B, we have 

1 ,.,6- 1 



a-D-2 



167T 



n a l b V [a (Sg b]c e) e 0l ... 



32tt 



-n 



l / fe V [a ^] 5g c c e 



a 1 ...a D _ 2 



(153) 



(154) 



Since 5g a b(, b = u5fi ab Z b + r5fi ab Y b , we have on B, 



2n a l b V [a (Sg b]c C 



n a l b (V a udnbcZ c - V b u5/j ac Z c + V a r5fi bc Y c - V b r5^ ac Y c ) 
l b 6fi bc Z c + n b 5/j bc Y c = , (155) 



as 5fi ab n a = = 5u Jab l a . Thus, the first term on the right side of SQ^\b is zero. To 
calculate the second term, we use, on B, 



2n a l b V [a £ b] 



n 



a l h [{Vaf)n b + {V a u)Z b + {V a r)Y b - {V b f)n a - {V b u)Z a - {V b r)Y a + 2/V [o n 6] ] 



n a Y a + l a Z a + n a V a f + -fl a V a {n b n b ) - fl b n a V a n b 
l a Z a + n a Y a + n a V a f 



(156) 



as n a V a r = n a n a = 0(r), and as V&(n a n a ) = V&(r 2 a) = 0(r), and n a V a n b = 0{r) since 
n a is tangent to affinely parameterized null geodesies on . Thus, we obtain 



(157) 



Thus, combining this with the previous result, we have shown that, with £ a = fn a + 
uZ a + rY a : 



[6Qi(g)-is6(g;6g)] 



1 

'8tt 



B 



f5# e - -(n a Y a + l a Z a + n a VJ)6e 



and, hence, 



S7T 



Wj,(g-6g,^g) = — / f5d e - ~(n a Y a + l a Z a + n a V a f)5e 



B 



+ c5H x 



(158) 



(159) 



It follows immediately from this formula that if 5g ab is a solution for which 5e\B = = 
and 6H X = 0, then W s (g;5g,^g) = 0. Conversely, if W^(g;8g,^g) = for all £ a of 
compact support, it follows immediately from ( 134"|) with E = that 8g ab satisfies the 
linearized constraints. It then follows immediately from (I159p that if W^(g; Sg,^g) = 
for all /, Z a , Y a , c, then 5e\ B = 6$\ B = 5H X = 0. □ 
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B Proof of lemmas [3] and [4] 

Lemma 3. Let Y = (N , Xq ) G Cq° © C^°. Then there exists a solution X in the space 
n fe (p 2 Wj +1 (E) © pWj(E)) C C°°(E) © C°°(£) to the following boundary value problem: 

CC*(X) = Y in E, (160) 

and X = (X, X a ) satisfies 

X a = Xrf , 5& = 5e = on 5 = dE. (161) 

Furthermore, if D > 5, the solution is unique, whereas for .D = 4, the solution is also 
unique unless t a is tangent to the generators of the horizon (i.e., the black hole is nonro- 
tating), in which case the solution is unique up to X a — > X a + ct a . 

Proof: The key ingredient in the proof is a Poincare-type inequality for certain tensor 
fields X = (X, X a ) G W 2 © W p l . Tensor fields X e W 2 ® W x p have, by the usual theorems 
about Sobolev spaces, a well defined restriction to B in the space W 3 ^ 2 (B) © W 1 ^ 2 {B) 1 
and the restriction map is continuous. Consequently, the space 

9£ ■= {X = (X, X a ) G W 2 p © W] I X a = N V a on B} (162) 

is a closed subspace of W 2 © W 7 ^. We first give the proof for D > 5 and then give the 
modifications to the proof for D = 4. 

D > 5: For X G we claim the following inequality in D > 5: 

C||^llvK p 2 ewi > ||£*M p (X)|| L 2 ffii 2 > e||X||vp2 eW n , M p = ^ ° j , (163) 

for non-zero constants c, C. The key inequality is the second one; the first inequal- 
ity is merely a straightforward consequence of the definition of £*, together with the 
asymptotic conditions fl27|) on the background spacetime and the fact that D a pD a p < 
C, p 2 D a D b pDbDbp < C, for other constants C . Combining the equality f)163p in the usual 
way with the Lax-Milgram theorem [T7] , we then infer the existence of X G 2£ satisfying 

= (L*M P X\£*M P ^) K + (Y\$) K for all $ G StT, (164) 

for any source Y G © in dimensions D > 5. Thus, taking Y = M p Y, and 
X = MpX, we get a distributional solution X to our equation (11601) in the space M P S£ . 

The solution is actually more regular, because, as is well-known [H [6], the operator 
CC* is an elliptic matrix operator of mixed type in the sense of Agmon-Douglis-Nirenberg; 
the leading derivative terms ("principal symbol") are in fact o~cc* — (X c* T(J c*j with X a = 
(X, X a ) and 




16. a c . [ Nc j = \^ pab _ ^_ habp 5{b ^ \\atc) (165) 

As a consequence, the standard techniques for elliptic operators (as adapted to systems of 
mixed type such as ours, see p. 210 of [28], and taking into account the weight factors of p) 
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then give that the solutions X are in fact in the Sobolev spaces X e Wp ®Wp 2 for all k, 
away from the boundary of S. Actually, this regularity extends to the boundary <9E = B 
by thm. 6.1 of [M], which is applicable in view of (I163p . By standard embedding 
theorems for Sobolev spaces, X, and also X, are in C°°(E) © C°°(E). It is then also 
straightforward to see that X solves the desired boundary value problem. First, because 
X e 3£ , we have N a = Nt] a , and because p is constant near B, the same holds for 
X a . Thus, we satisfy the first of the required boundary conditions. Next, because X is 
already known to be a weak solution of the equation (I160p . and because it is smooth, 
it must satisfy the equation in the classical sense. Furthermore, performing a partial 



integration in eq. ( I164p , with X 



M p X,Y 



M X F, and smooth if) compactly supported 



in E, then shows that be boundary terms must vanish. The boundary terms are obtained 
from ( I78p . Because if) has compact support in E, only a boundary term from B arises, 



which, with (5h, dp) :- 
In view of eq. 



C*MpX, is given by = j B [5Q${g) — i$0(8g)], where if) = M p if>. 
1158ft this implies that = = 5e\s- These are precisely the last two 
remaining boundary conditions in (1161ft . 

Thus, what remains is to demonstrate the weighted Poincare- inequality (1163ft . We 
first consider smooth X with compact support in a compact set K C E, where p may 
be assumed to be 1. In other words, such X are zero in the asymptotic region, but not 
necessarily on B. For such fields, we can find a constant c > such that 



/ h^(D b N a N b -D b N b N a )r 1 a 

JB 



(166) 



For a proof, see lemma 2.6 of [5J. Now we assume that N a = Ni] a on B. This means 
that, near B, we can write iV a = Nr] a + xY a , with x = dist# the geodesic distance to B, 
and Y a some smooth tensor field. Then the boundary term becomes 



/ h? (D b N a N b - D b N b N a )r] a 

JB 



hi (ND b Nr] b + N 2 r] a r] b D b r] a + NY a r] a - ND b Nr] b - N' 2 D a r] a - NY a r) a ) 



B 



N 2 H hh = 



(167) 



B 



where H is the trace of the extrinsic curvature of B inside E. To simplify the calculation, 
we have assumed that r] a is extended away from B such that 7] a D a 7] b = 0, and we used 
D a x = T] a . Because B is the bifurcation surface, H = 0, and the boundary term vanishes. 
Hence, (I166I) holds without the boundary term if X satisfies iV a = Ni] a on B, i.e. if 

x e % n (C£°(K) © C%°(K)). _ 

Consider next an open region O C E such that O C K. It is shown combining Prop. 7.2 
and Lemma 4.2 of @, that for all X e C§°(£ \ O) © C °°(E \ O), we have 



c\\X\ 



W 2 OW 



x < ||£*M p (X)|| L2eL2 



(168) 



for some constant c > 0. We now combine this and eq. (11661) (without boundary term) 
to obtain the desired inequality (I163p . Suppose, to obtain a contradiction, that there 



37 Note that the arguments in this paper are entirely local near B, so it does not matter that our domain 
is not bounded. 
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is a sequence X n G JT, such that ll^nllw 2 ©^ = 1 3 but ||£*M p (X n ) ||l 2 ©l 2 — ► 0. Let 
Xi + X2 = 1 be a partition of unity such that suppxi C K, suppx2 CS\0, Then we 
estimate 



ll^llwlffiWp 1 < IIXl^nlllV^M/i + ||X2X"n||w2©W"i 

< CH^^Xi^n)^ 2 ®^ + C||xiX n ||v[/i©w + C\\C* M p (x2X n )\\ L 2 eL 2 

< C||Xl£*(X„) II L 2 ©L 2 + C||[£*,Xl]X n || L 2 ei 2 + CHxiXnUtyi^O + 

C\\ X 2C*M p (X n )\\ L2(SL2 + C\\[C*, X 2]M p X n \\ L2(SL2 

< C||£*M p (X n )|| L 2( S ) eL 2( 2 ) + C||X n ||w-i(o)0W°(o) (169) 

with possibly new constants in each line. In the last step we used that the commutator 
[£*, if)] with a smooth compactly supported function if) decreases the order of each entry of 
the matrix operator C* by one unit (unless the order of the entry is already = 0), so that 
[£*,ip] : W 2 © W 1 -> W 1 © W° is bounded. Now, by assumption, ||£*M p (X„)|| L 2 eL 2 -> 
for the first term on the right side. On the other hand, since X n is by assumption bounded 
in W 2 (0) © W l (0), it follows from the Rellich-Kondrachov compactness theorem (see 
e.g. [E]) that X n (or a subsequence thereof) is Cauchy in W x {0) © W°(0). Hence, the 
above inequality shows that a subsequence of X n is Cauchy in 3£ C W 2 (E) © (£), 
hence convergent with limit X G Since the norm of X n in this space is = 1, X ^ 0, 
and by the continuity of C*M p : -»■ L 2 © L 2 , we learn that ^ := M p X G p 2 ^ 2 © pW 7 "^ 
satisfies £*(£) = 0. So ^ a must be equal, almost everywhere, to a non-trivial Killing 
vector field. However, a nontrivial Killing field in an asymptotically flat spacetime must 
approach an infinitesimal Poincare transformation near infinity. Consequently, in D > 5 
spacetime dimensions, there are no non-zero Killing for which M~ x £ G L 2 © L 2 , hence a 
contradiction. 

D = 4: In this dimension, the last statement is not true for the timelike Killing field £ a = t a 
if the background is such that t a = (N, N a ) is tangent to the generators of the horizon: 
Indeed, M~ x t then is in L 2 © L 2 , and it satisfies the boundary condition N a = Nr) a . 
Consequently, the Poincare inequality (I163p is not true in D = 4 e.g. for a Schwarzschild 
background. However, in that case we can modify our argument as follows. We now 
supplement the definition of the space 3^ eq. f)162p by the condition that elements of 5£ 
should be L 2 ©L 2 orthogonal to M~ x t. Then, in this space, the Poincare inequality holds, 
and the rest of the argument goes through, giving a unique solution to C*C(X) = Y 
in Mp^fc with the desired boundary conditions. The general solution is then simply 
X a + ct a . □ 



Lemma 4. Let j a be a smooth vector field on E such that j a G p x W p 1 for any k > 0. 
Then there exists a unique, solution X a G pW p to the boundary value problem 

D a (D {a X h) - j^-jh ab D c X c ) = / (170) 

and 

_ v a V b )X a \ B = , r) a D a ( V b X b )\ B = . (171) 
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Furthermore \\p X a \\ W k < C\\pj a \\ w k-i. 

The same statement is true for the boundary value problem for X a on the slice S = 
S x F in the black brane spacetime with spatial metric h ab = h a b + ^2(d,Zi) a (dzi)b, with 
D a replaced by D a , and with the coefficient in the trace term in the operator replaced by 
l/(D + p-l). 

Proof: Define (PX) ab = D^ a X b ) — jy—^h a bD c X c . The key ingredient of the proof of is the 
weighted Poincare inequality 

c / p~ 2 X a X a hh < [ (PX) ab (PX) ab hh (172) 
Jt Jt 

for some c > 0. This inequality holds for any smooth X a satisfying the boundary condi- 
tions stated in lemma HJ In order to prove this inequality, we consider the identity 

(pxy b (px) ab = 

1 1 D — 3 ( 173) 

- (D a X b )D a X b + D a (X^D b X^) - - R ab (h)X a X b + ^ _ - (D C X C ) 2 , 1 } 

which holds for any vector field. Now let O = {p < C} for some large C, and assume 
that X a has compact support in E \ O. Integrating the above identity then gives 

^J^D a x b )D a x b hh < J^(pxy b (px) ab hh + 1 J R ab (h)x a x b 

< [ (PX) ab (PX) ab hh + isupp 2 (R ab (h)R ab (h)fz [ P - 2 x a x a hh 



< / (PX) ab (PX) ab hh + e j^p- 2 X a X a hi 



(174) 



for any given e > 0, provided C is chosen sufficiently large (here we use that R a b(h) = 
0(p~( D ~ 1 - ) )). Now from prop. C.5 of |6J, we also have for sufficiently large C 

c [ p- 2 X a X a hh < [ (D a X b )D a X b hh , (175) 

JT, JT 

and some c > 0. It follows that inequality (I172p holds if X a is supported outside a 
sufficiently large O. Consider next a compact set K such that K D O, and let X a now 
be smooth and supported in K. Integrating now our identity ( 1 1 73 [) over K, we obtain 



c\\X\\ m < \\PX\\ L 2 + \\X\\ L 2 + 



h*{D b X a X b -D b X b X a )rf 



B 



(176) 



for some c > 0. Let now = {X a e pW p \ = {h ab - p a p b )X b \ B }. Then the same 
calculation as in (I167P shows that the boundary term vanishes for X a G 3C . Thus, for 
X a £ 1" we have eq. (I172p if the support of X a is outside O, whereas we have (11761) 
without boundary term if the support of X a is in K. By considering now a partition 
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of unity for the covering O U K = E, and arguing precisely as above around (I169p . we 
conclude that the Poincare inequality (11721) holds for any X a G E£ , unless there is an 
X a G 3£ such that (PX) a b = 0. By the definition of P, such a X a would be a conformal 
Killing vector field for h a b, and would hence have to approach a conformal Killing vector 
field on R-° _1 in the asymptotic region. In particular, the slowest possible fall-off of X a 
would be that of a translation. However, by the definition of SC , we would also have to 
have p~ 1 X a G L 2 , which is impossible. This establishes that the Poincare inequality holds 
for any X a G SC. 

Combining the Poincare inequality in the usual way with the Lax-Milgram theorem 
on 5C (see lemma [3] for a similar argument), we get a weak solution X a G S£ to our 
boundary value problem in the weak sense, i.e. (Pip\PX) = —(ip\j) for all ip a £ S£ ■ This 
also gives H/T^Jf < C 1 1 PJ 1 1 - The usual bootstrap and scaling arguments then show 
< Cllpj'll^fe-i for any k, so the weak solution is a strong solution satisfying 

the Neumann condition i] a D a (r] b Xb)\B = 0, whereas the remaining boundary conditions 
are already built into SC . This completes the proof of the lemma in the asymptotically 
flat case. The asymptotically KK case is exactly the same. □ 



C Strengthened version of property (3) of Proposi- 
tion E 

In this Appendix, we strengthen property (3) of proposition [5] with the following lemma: 

Lemma 5. Let (Sh, Sp) be a perturbation in the dense subspace V D (Wp +1 © W^) for 
sufficiently large k. Then for any S > 0, we can find a new perturbation (Sh', Sp') in the 
same subspace having the property that it vanishes in a neighborhood of spatial infinity, 
and is such that \\(6h,6p) - (Sh', 5p') \\ w *>+^ W k < S. In fact, if (Sh,Sp) G V H (C°°(S) © 

C°°(S)), then (Sh', Sp') can be chosen in Vn (C£°(E) ffiC °°(E)), i.e. smooth and vanishing 
near spatial infinity. 

Proof: The proof of this lemma follows from the (much more general) techniques and 
results in [HE], and also [6]. We consider a sufficiently large R > and the corresponding 
annular domain Q = {p G £ | R < p(p) < 2R}, with p > as usual a function which 
is equal to the radial distance in the asymptotic region. Furthermore, we introduce a 
function x : Q — > R which in a small neighborhood of dQ is given by x = distgn, and 
satisfies x > in the interior. A new norm for functions or tensor fields u G C^°(fl) is 
then introduced by 



u 



w 



• • • D an) u)D^ ■ ■ ■ D a ^u (p x 2 ) 2n e- 2 ^ J , (177) 



and W!? 2 t /_ is defined to be the closure of such functions or tensor fields in this norm. 

0,px ,e L ' x 

Now let x be a smooth function which is equal to 1 in the interior of E, which interpolates 
between 1 and strictly inside Q, and which is for p > 2R. We can choose this function so 
that (D ai ...D an x)D ai ...D an x < CR~ 2n . Our first try for the compactly supported initial 
data is (\Sh,\Sp\ Then for each S > 0, k G N we can find a sufficiently large R so that 
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\\(5h,5p)-(x5h, X Sp)\\ 

w$ +l ®w k < Of course, (xSh, xSp) does not satisfy the linearized 
constraints; £{xSh,x$p) =: Y is a non-zero element Y G W^ 2fi _ 1/x (Q) © W^ 2 ^_ 1/x (fi) 
(with the usual decomposition of Y a into lapse and shift understood). To satisfy the 
linearized constraints, we wish to correct {xSh, x$p) by a suitable element of the form 
—e~ 2 ' x M 2 £*(X), with M x = diag(x 2 ,x) as above, i.e. we wish to set 

(Sti, Sp') := (x6h, X 5p) - e- 2 ' x MlC*{X) (178) 

The factor e~ 2 ^ x is inserted in order that we can continue the second term sufficiently 
smoothly by to all of S. Such a perturbation is then clearly compactly supported, and, 
to satisfy the linearized constraints, X would have to satisfy 

e V*£e-*/*M*£*(X) = e 2/x Y . (179) 

The operator e 2 ^ x £e~ 2 ^ x M 2 £* on the left hand side is shown in [6] to be a bounded map 

where £ = span(t a , ipf, ...,if)ff) is the span of the Killing fields, and _L means orthogonal 
complement in L 2 (Q,e~ 2 ' x h^). Furthermore, the operator has a bounded inverse. It 
follows from SM, SJa = for the perturbation (Sh, dp) G V [cf. Proposition [5] that 

<£i e2/ * y >L W Vv*) = (ZWh, X Sp))i>W = (180) 
(£*(0\Y)m m mx) - SH^Sh, 5p) = + 

for any £ G t. Hence e 2 l x Y is in Furthermore, for sufficiently large R, the x2 e _ 1/x (fi)© 

2 e _ 1/a .(fi)-norm is arbitrarily small. So we can invert ( I179P with X G + Jj e _ 1/x © 
W^ e _ 1/X having small norm, and hence with -e~ 2/x M 2 £*(X) G e~ 2/x M 2 (W^ 2 e '_ 1/a: (fi)© 
Wq x2 e _ 1/x (Q)) having small norm in W / ^ +1 (E)©H / ^ : (E), where we extend the tensor fields 
by outside fl Consequently, for sufficiently large R, the Wp +1 (T>) © VF*(£)-norm of 

(jit) ~ {tp) = (& ~ % H p) " ^ G ® ^ ( 181 ) 

is as small as we like. By construction (Sh', Sp') satisfies the linearized constraints, and 
has SM = SP = S3 = and S^Ib = Se\s = 0. Thus, it is in the space V by proposition |5j 
As shown in prop. 5.7 of [B], if (Sh,Sp) G V is even smooth, then — e~ 2 ' x M 2 £*(X) is 
a smooth pair of tensor fields on Q, which can be continued by outside of Q and the 
resulting tensor fields (Sh',Sp') are in Cq°(£). □ 
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